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ABSTRACT 

The  existence  of  multi-mode  propagation  can  have  two  deleterious 
effects  on  an  OTH  radar  tracker;  it  can  lead  to  ambiguity  in  the  interpre- 
tation of  the  range  of  detected  targets,  and  it  can  lead  to  difficulties  in 
determining  their  number.  Because  the  identity  of  a mode  cannot  be  uniquely 
determined  from  a single  detection  sample,  any  testing  to  determine  mode 
identity  must  be  performed  at  or  beyond  the  tracker. 

Three  methods  are  considered  for  identifying  propagation  modes; 
one  based  on  testing  for  the  occurrences  of  echoes  in  regions  or  combina- 
tions of  regions  of  the  radar-observable  space  accessible  to  only  a single 
propagation  mode,  another  based  on  the  amplitude  fluctuations  observed  on 
echoes,  and  a third  in  which  mode  separations  are  tracked  and  taken  into 
account . 

Predictive  filters  are  given  general  treatment.  The  Kalman 
filter  tracking  algorithm  is  discussed  with  the  modifications  required 
to  implement  the  first  and  third  (above)  mode-identification  techniques. 

The  use  of  signal  strength  fluctuations  ss  a mode-identification  criterion 
is  examined  and  dismissed  as  being  insufficiently  predictable  to  provide 
a reliable  means  for  identifying  propagation  modes. 
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(Uj  The  investigations  reported  herein  were  conducted  during  the  period 

February  through  December  1974  under  the  sponsorship  of  the  Office  of  Naval 
Research,  Code  464,  under  Contract  N00014-74-C-0106. 

1 . 1 The  Role  of  Mode  Identification  in  an  OTHR  System 

To  permit  OTH  radar  data  to  be  interpreted  and  displayed 
in  geographic  coordinates,  a coordinate  conversion  process  is  required. 

Because  the  conversion  process  must  be  done  differently  from  one  mode  to 
another,  it  is  necessary  to  correctly  identify  the  propagation  mode  by 
which  each  target,  is  detected.  In  some  cases,  this  is  no  problem  because 
by  properly  choosing  frequency  and  other  operating  parameters  of  the  radar, 
it  is  possible  to  assure  that  targets  are  illuminated  via  a unique  mode. 
However,  all  too  often  such  a solution  is  precluded,  either  by  the  impossibil- 
ity of  managing  propagation  so  as  to  isolate  a single  mode,  or  simply  because 
conflicting  considerations  make  such  a choice  of  operating  parameters  inadvis- 
able. 

In  particular,  propagation  management  is  already  greatly  constrained 
by  the  need  to  provide  adequate  illumination  in  the  surveillance  area,  and 
it  would  be  undesirable  to  compromise  the  amount  of  energy  placed  on  target 
for  other  considerations.  The  methods  to  be  discussed  later  are  intended 
to  avoid, or  at  least  to  minimize,  the  need  for  such  a compromise. 

There  is  no  practicable  way  in  which  the  mode  of  propagation  can 
be  determined  from  a single  pulse  or  single  doppler  dwell  observation  of 
a target.  Although  in  principle  such  a determination  could  be  made  by  measur- 
ing the  vertical  angle  of  arrival  of  the  echoed  radiation,  in  practice  this 
would  require  greater  antenna  resolutions  than  are  economically  feasible. 

Most  of  the  methdds  to  be  considered  here  involve  multiple  observation  of 
the  target--that  is,  observations  made  during  different  pulses  or  doppler 
dwells  of  the  radar.  Because  such  multiple  observations  of  a target  are  not 
associated  until  the  "tracker"  stage  of  the  radar  data  processing,  the  methods 
to  be  discussed  here  are,  of  necessity,  applied  at  or  beyond  the  tracker.  One 
technique  that  has  been  considered  can  be  applied  within  a single  frame  of 
data  but  can  be  applied  to  tracked  data  as  well.  Therefore,  the  discussion 
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will  hereafter  assume  that  tracked  data  are  to  be  used.  In  such  a mode  of 
operation,  the  tracker  operates  on  data  in  terms  of  the  radar  cbservables-- 
time  delay  of  the  echo,  doppler  shift  of  the  echo,  and  angle  of  the  echoed 
radiation  relative  to  a reference.  The  only  requirements  imposed  on 

the  tracker  are  that  it  be  tolerant  to  any  discontinuities  of  the  track  and  of 
apparent  accelerations  of  the  targets  produced  by  the  propagation  character- 
istics of  the  ionosphere. 

1 • 2 The  Problems  of  Ambiguous  and  Multiple  Returns 

For  the  purpose  of  the  present  discussion,  it  is  useful  to  identify 
two  parts  of  the  problem  caused  by  the  radarist's  uncertainty  regarding  the  mode 
of  propagation.  In  order  to  correctly  interpret  the  data  from  an  OTH  radar,  it 
is  necessary  to  perform  a one-to-one  mapping  of  the  target's  position  in  radar 
observable  space  to  geographical  coordinates.  However,  under  multiple  mode 
conditions,  this  mapping  process  is  complicated  in  two  ways.  First,  the 
detected  target  may  occur  at  a point  in  the  radar  observable  space  that  can 
be  mapped  to  any  of  two  or  more  points  in  the  true  geographic  space,  depending 
upon  what  mode  is  assumed  to  have  supported  the  propagation.  Until  the  re- 
sponsible mode  can  be  identified,  it  is  necessary  for  the  radarist  to  carry 
along  a set  of  multiple  hypotheses  regarding  the  position  of  the  detected 
target.  In  effect,  the  data  must  undergo  a one-to-many  mapping  from  the 
radar- observables  space  to  geographic  coordinates.  We  here  speak  of  this  as 
the  mode  ambiguity  problem,  the  effects  of  which  can  be  stated: 

1.  The  true  number  of  resolved  targets  is  readily  apparent. 

2.  But  the  true  position  of  the  target  is  one  of  perhaps 
several  discretely  determined  positions. 

The  second  kind  of  problem  occurs  if  the  target  happens  to  be  at  a 
point  in  the  space  of  geographical  coordinates  at  which  it  is  illuminated 
via  more  than  one  propagation  mode.  The  effect  of  this  will  be  that  the 
target  will  manifest  itself  at  several  points  in  the  radar- observable  space. 

A distinct  return  from  the  target  can  occur  for  each  propagation  mode. 

Additional  returns  will  also  arise  corresponding  to  outward  propagation  via 
one  mode  and  return  propagation  via  another.  However,  mixed-mode  returns 
will  occur  with  the  same  delay  regardless  of  the  order  in  which  any  given  pair 
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of  modes  occur- -whether,  for  example,  an  E mode  occurs  in  the  outward 
leg  and  an  F mode  on  the  return  or  vice  versa.  In  view  of  the  above,  we 
see  that  the  number  of  returns  from  a target  in  a region  of  N modes  of 
illumination  will  be  equal  to  the  number  of  ways  N elements  can  be  drawn  twice, 
with  replacement,  to  form  unordered  pairs.  Thus,  the  number  of  returns,  n,  is: 

N 

n = £ i 

i-1 

Unless  the  radar  has  resolution  on  the  order  of  1 millisecond, 
the  ordinary  and  extraordinary  modes  of  propagation  will  remain  unresolved 
and  manifest  themselves  only  as  a polarization  rotation  and  associated 
fading.  Propagation  via  the  F layer  will  occur  in  two  modes--the  high  ray 
and  the  low  ray.  In  addition,  there  exists  the  possibility  «f  propagation 
via  the  E layer  and  FI  layers  when  they  exist.  Depending  on  the  height  and 
electron  densities  of  the  layers,  various  combinations  of  these  may  propa- 
gate radar  illumination  to  the  region  of  surveillance  and  give  rise  to  the 
multimoding  effects.  Note  that  the  effect  of  multimoding  differs  from 
that  of  mode  ambiguity  in  that: 

1.  The  true  number  of  targets  cannot  be  determined  until 
an  association  is  made  among  the  multiple  echoes  pro- 
duced by  each  individual  target. 

2.  Knowledge  of  the  true  position  of  the  target  is  not 
precluded  by  multiple  moding.  Only  if  one  or  more  of 
the  returns  lies  in  ambiguous  regions  of  the  radar- 
observable  space  will  there  be  range  ambiguity  of  the 
kind  discussed  earlier. 

The  practical  effect  of  multimoding  is  that  it  may  produce  erroneous  estimates 
of  the  number  of  targets  detected  in  the  surveillance  region.  A correct 
coordinate  conversion  from  radar  observable  space  to  geographic  coordinates 
would  lead  to  coincident  points,  readily  interpreted  as  being  a single 
target.  However,  all  too  often,  errors  in  the  coordinate  conversion  process 
will  lead  to  imperfect  coincidence  and  difficulty  in  interpretation. 

A second  practical  implication  of  the  multiple  mode  effect  is 
that  the  multiple  returns  could  be  combined  non-coherently.  or  perhaps  even 
coherently,  to  improve  the  detection  sensitivity  of  the  radar.  However,  to 
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do  this,  it  would  be  necessary  to  convert  from  radar  observable  coordinates 
to  geographical  coordinates  prior  to  tracking,  and  various  hypothetical  tracks 
would  have  to  be  carried  for  radar  returns  occurring  in  ambiguous  regions  of 
the  radar  observable  space.  Consequently,  although  better  performance  could 
be  expected  with  such  an  arrangement,  it  would  be  at  the  expense  of  more 
powerful  data  processing.  The  coordinate  conversion  processor  would  have 
to  be  powerful  enough  to  perform  a conversion  for  each  ambiguity  of  each 
detection,  and  the  tracker  would  have  to  be  powerful  enough  to  carry  the 
additional  hypothetical  track.,  that  arise  from  the  mode  ambiguity. 

1.3  Mode  Identification  and  ECM 

Many  of  the  considerations  involved  with  mode  identification  also 
have  applicability  to  the  problem  of  making  the  radar  immune  to  deceptive 
forms  of  electronic  countermeasures.  In  order  to  produce  false  targets  in 
an  OTH  radar,  it  is  necessary  to  transmit  a false  pulse  delayed  with  re- 
spect to  the  time  of  transmission  of  the  radar  illumination  pulse  by  an  amount 
appropriate  to  the  range  at  which  the  false  target  is  to  appear.  However, 
a sophisticated  radarist  employing  mode  identification  techniques  would  not 
depend  only  upon  the  delay  of  the  radar  signal  as  an  indication  of  the 
apparent  target  range.  He  would  also  seek  other  cues,  such  as  the  maximum 
frequency  illuminating  the  target  via  a particular  mode,  and  the  degree  of 
coincidence  among  positions  of  a target  as  deduced  from  multiple  modes  of 
coverage.  The  need  for  sufficient  information  to  simulate  all  of  these 
effects  would  make  the  task  of  the  radarist 's  ECM  opponent  much  more 
difficult,  and  were  proper  measures  taken  on  the  part  of  the  radarist,  vulner- 
ability to  spoofing  ECM  could  nearly  be  precluded. 

Consider  a radar  operating  with  a nonrepet i t ive,  unpredictable  wave- 
form- -such  a radar  is  susceptible  only  to  false  targets  at  delays  greater 
than  that  of  the  spoofing  transmitter.  (For  the  moment,  we  assume  that  we  are 
dealing  only  with  a single  propagation  mode.)  However,  if  the  l-adarist 
increases  frequency  until  the  point  is  reached  at  which  the  spoofing  transmitter 
(and  receiver)  is  in  the  skip  zone,  while  the  apparent  position  of  the  target 
is  still  being  illuminated,  the  spoofing  can  no  longer  be  carried  on  because 
the  spoofer  becomes  decoupled  from  the  radar.  The  ability  to  keep  track  of 
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tests of  this  kind  is  a large  part  of  the  mode  identification  process,  so 
the  ability  to  recognize  spoofing  ECM  is  an  additional  benefit  to  be  gained 
by  employing  the  techniques  discussed  here. 

I*4  Fixed  Frequency  Versus  Multiple  Frequency  Operation 

Many  of  the  techniques  here  depend  on  the  radar  having  some 
degree  of  frequency  agility.  For  our  purpose,  it  will  usually  be  satis- 
factory for  the  radar  to  be  able  to  change  frequency  from  one  doppler  dwell 
to  the  next.  The  changing  of  frequency  can  either  be  carried  out  in  a 
predetermined  pattern  (e.g.,  by  assigning  frequencies  by  means  of  a pseudo- 
random sequence)  or  can  be  made  adaptive  to  the  set  of  echoes  of  the  moment. 
The  latter  procedure  would  be  the  preferred  one  if  the  frequencies  required 
to  identify  propagation  mode  were  unsatisfactory  for  general  illumination 
of  the  surveillance  area.  However,  such  a procedure  would  require  that  the 
tracker-mode  identifier  be  given  the  capability  of  selecting  the  radar’s 
operating  frequency. 
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If  the  radar  is  constrained  to  fixed- frequency  operation,  the 
ability  of  the  radar  to  carry  out  tests  to  determine  propagation  modes  is 
considerably  limited.  Tests  involving  determination  of  the  maximum 
observable  frequency  to  a target  then  reduce  to  observations  of  the  radar 
delay  at  which  the  target  enters  the  mode's  skip  zone.  Because  this  kind 
of  test  is  not  likely  to  be  generally  satisfactory,  consideration  was  given  to 
a second  set  of  mode  identification  tests  based  on  the  time  variation  of  the 
amplitude  of  the  radar  echo  from  the  target.  A discussion  of  this  possibility 
is  to  be  found  in  Section  4. 
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2.  THE  PROBLEMS  OF  MODE  AMBIGUITY  AND  MULTIPLICITY 


2. 1 Prevalence  of  Multimode  Conditions 

As  mentioned  in  Section  1,  it  is  often  possible  to  take  steps 
in  the  management  of  the  OTH  radar  to  reduce  the  number  of  modes  with 
which  the  radar  must  contend.  A good  example  of  this  is  the  practice 
of  restricting  the  radar's  operation  to  radiation  propagated  via  the 
F modes.  This  is  done  by  choosing  a sufficiently  b:.gh  operating  fre- 
quency so  that  the  waves  penetrate  the  E layer  except  at  very  low 
radiation  angles.  At  these  low  angles  the  combined  effects  of  reduced 
antenna  gain  and  D- layer  absorption  largely  reduce  the  radar's  response 
to  clutter  and  targets  at  ranges  of  interest. 

Unfortunately,  there  are  many  circumstances  in  which  this  is  not 
possible;  for  example,  when  the  critical  frequencies  of  two  layers  are 
not  sufficiently  different.  F.ven  when  the  effects  of  other  ionospheric 
layers  can  be  avoided,  it  is  still  necessary  to  contend  with  the  various 
forms  of  multiraoding  within  the  layer.  The  occurrence  and  character- 
istics of  the  various  forms  of  multiraoding  are  discussed  next  for 
background. 

2.1.1  High  Ray- Low  Ray  Propagation 

Each  ionospheric  layer  supports  radio  propagation  between  a 
pair  of  points  via  several  modes.  :-,ome  of  the  multiplicity  of  modes  are 
the  result  of  the  earth's  magnetic  field  on  propagation.  However, 
ignoring  for  the  moment  the  mode  splitting  attributable  to  magneto- 
ionic action,  we  are  left  with  two  propagation  modes  that  are  supported 
by  each  layer:  the  low  angle  ray  and  the  high  angle  or  Pedersen  ray. 

It  is  the  low  ray  that  is  normally  counted  on  to  provide  the  propagation 
from  radar  to  target,  being  more  stable  and  not  subject  to  the  defocusing 


action  to  which  the  Pedersen  ray  is 
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The  high  and  low  ray  modes  of  propagation  can  produce  three 
signals  from  each  target:  one  trace  produced  by  echoes  propagated  by 

the  low  ray  mode,  a second  produced  by  high  ray  propagation,  and  a 
third, at  intermediate  delay,  which  is  produced  by  echoes  that  have  propa- 
gated by  high  ray  in  one  direction  and  low  ray  in  the  other. 

In  this  case,  not  only  is  there  a multiplicity  of  modes  but  also 
an  ambiguity  because  each  trace  lies  in  a region  of  the  frequency- delay 
plane  that  is  accessible  to  the  other  two  signal  paths.  Although  the 
association  of  the  traces  is  readily  apparent  on  an  F scan, where  they 
join  at  the  layer's  skip -zone  caustic,  their  association  would  not 
be  as  readily  apparent  in  a single  radar  observation  or  in  a series  of 
observations  at  a single  frequency. 

For  strong  echoes,  such  as  those  from  a repeater,  two  or  more  of 
these  modes  can  often  be  seen  to  extend  several  megahertz  below  the 
junction  frequency.  However,  the  high  ray  is  subjected  to  an  appreciable 
defocus ing  effect,  causing  signals  propagatdd  by  that  mode  to  be  sub- 
stantially weaker  than  those  propagated  via  the  low  ray,  particularly 
as  the  operating  frequency  departs  more  and  more  from  the  junction  fre- 
quency for  the  target's  range. 

Tornatore  (1972)  of  ITT-EPL  has  measured  the  relative  strength  of 
signals  propagated  one  way  via  the  high  and  low  rays  and  has  found  that 
the  low  ray  signal  strength  exceeded  that  of  the  high  ray  by  median 
ratios  of  3,  6 and  10  db  as  measured  1,  2,  and  3 MHz  below  the  junction 
frequency,  respectively.  These  results  can  be  applied  directly  to  the 
two-way  radar  case  to  give  the  relative  amplitude  of  the  low  and  mixed- 
ray  modes  but  must  be  doubled  to  give  the  ratio  of  low-to-high  ray 
radar  signal  strength. 


Figure  2-1  shows  an  oblique  ionogram  indicating  radar  propagation 
via  the  F la>„r  between  Whitehouse,  Virginia,  and  Lincoln,  Nebraska.  The 
regions  of  high,  low  and  mixed-mode  ambiguity  are  indicated  in  the  lower 
tracing. 
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2.1.2  Multiple  Layer  Multimoding 

Under  daytime  conditions  the  ionosphere  frequently  comprises 
two  or  more  layers  capable  of  supporting  propagation;  and  even  under 
nighttime  conditions,  during  spring  and  summer  months,  a sporadic  E 
layer  frequently  accompanies  the  F layer  and  produces  multimoding. 

Depending  on  the  distribution  of  electron  densities  in  the  layers, 
either  the  F mode  or  an  underlying  mode  is  capable  of  supporting  propa- 
gation to  a given  distance  at  higher  frequencies.  The  anomalously  high 
F layer  critical  frequency  occurring  on  winter  days  are  typical  of  those 
cases  where  the  F layer  supports  the  higher  frequency  of  propagation, 
permitting  undesired  modes  to  be  avoided  by  selecting  a suitably  high 
operating  frequency.  On  the  other  hand,  during  summer  days,  the  opposite 
condition  exists  and  the  E layer  often  supports  the  higher  frequency  of 
propagation  to  any  range.  Figure  2-2  shows  the  computed  traces  of  leading- 
edge-focused  ground  clutter  superimposed  on  oblique  ionograms  to  two 
jes.  These  traces  were  computed  for  a model  ionosphere  appropriate 
j summer  daytime  propagation  southward  from  the  ONR,  Whitehouse,  Virginia, 
site.  Note  that  the  E layer  is  shown  to  support  propagation  to  radar 
delays  of  8 to  16  ms  at  frequencies  some  4 MHz  higher  than  does  the 
F layer.  However,  the  radiation  angles  required  to  achieve  the  E layer 
propagation  are  quite  low,  being  about  seven  degrees  for  1000  km 
propagation  (=8  ms)  and  going  to  zero  degrees  before  the  16  ms  delay 
can  be  achieved.  The  reduced  gain  of  practical  antennas  at  such  low 
radiation  angles  usually  prevents  E-mode  propagation  except  to  relatively 
short  ranges.  Consequently,  in  these  cases  the  E layer  does  not  directly 
contribute  to  the  mode  ambiguity  and  multimoding  problems.  However,  it 
does  have  a serious  indirect  effect — it  obscures  the  higher  layers 
over  a broad  range  of  frequencies.  This  effect  is  shown  in  Figure  2-2 
by  the  abruptly  truncated  low  ray  traces.  The  result  is  that  the  radar 
must  be  operated  at  a frequency  near  the  MUF  of  the  F layers.  This 
leads  not  only  to  the  necessity  to  cope  with  the  high  ray  propagation  but 
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Traces  of  Lead ing-Hdge- Focused  Ground  Clutter  and  Discrete 
Targets  at  Ranges  of  2554  and  1063  KM,  Computed  from  a Model 
Ionosphere  Typical  of  a Summer  Daytime.  Propagation  southeast 
from  Whi rehouse,  Virginia,  is  assumed. 


also  with  a multimoding  problem  involving  FI  and  F2  layer  propagation. 
Figure  2-3  is  a skywave  backscatter  record  made  with  a large  cross 
section  repeater  in  the  field  of  view.  Note  that  the  repeater  trace, 
which  has  been  produced  by  F layer  propagation,  does  not  have  a low- 
ray  component  extending  to  low  frequencies.  Rather,  it  is  limited  to 
the  frequencies  near  the  MOF  where  the  high  ray  and  mixed  ray  produce 
raultimoding. 

2.1.3  Sporadic  E .'lodes 

During  the  sp*ing  and  summer  (near  summer  solstice)  months,  an 
additional  ionospheric  component  compounds  the  multimoding  problem 
further.  The  E layer  is  subject  to  localized  enhancements  leading  to 
critical  frequencies  that  can  be  in  excess  of  7 MHz.  Beside;  ntensi- 
fying  the  obscuration  problem  posed  by  the  normal  E layer  during  the 
daytime  hours,  sporadic  E can  also  introduce  similar  problems  at  night. 
The  probability  of  its  occurrences  with  critical  frequencies  in 
excess  of  5 ffiz  is  on  the  order  of  30%.  Figure  2-4  shows  a skywave 
backscatter  record  in  which  sporadic  E is  strongly  evident  at  ranges 
out  to  about  10  msec.  In  these  cases,  propagation  via  the  sporadic  E 
layer,  although  limited  and  uncertain,  can  be  managed  to  avoid  multi- 
moding  effects.  The  same  cannot  be  said  of  F-mode  propagation.  Simple 
sporadic  E bkk  s will  not  pose  a problem  at  radar  delays  in  excess  of 
10  msec.,  but  multiple  hop  modes  are  a possibility. 

The  likelihood  of  two-hop  Es  modes  is  reduced  by  the  requirement 
that  sporadic  E must  occur  at  two  independent  reflection  points.  How- 
ever, if  the  operating  frequency  is  chosen  too  far  below  the  F layer 
MOF  for  the  range  of  interest,  the  possibility  of  two-hop  modes  involving 
the  E and  F layer  arises.  The  probability  of  these  occurring  is 
approximately  the  same  as  for  the  one-hop  E mode,  since  only  one  E -layer 
reflection  is  required. 
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figure  2-5.  Skywave  Iiacksc.it  ter  f-Scan  Showing  He  ho  Supported  by  l-l.i>or 
Propagation  Which  is  Confined  to  frequencies  \ear  the  MIH 
Because  of  liiuicr lying  Obscuring  layers 


2.1.4  Multihop  Modes 


Because  multihop  modes  Tequire  ionospheric  reflections  at  angles 
more  nearly  vertical  than  do  single  hop  modes,  there  is  usually  an  in- 
terval of  frequency  at  which  the  single  hop  mode  is  supported  but  which 
penetrates  the  layer  at  the  multihop  angle.  Figure  2-5  shows  an  ex- 
ample of  a one-hop  mode,  a two-hop  mode,  and  the  intermediate  mixture 
mode.  Note  that  the  mixture  mode  and  the  two-hop  mode  both  reach  an 
upper  frequency  limit  of  about  8 MHz  lower  in  frequency  than  does 
the  one-hop.  This  interval  would  provide  ample  scope  for  avoiding 
multimoding  problems. 

2.1.5  Unresolved  Modes 

In  addition  to  the  multiple  moding  effects  described  above,  there 
are  other  forms  of  multimoding.  The  splitting  of  propagated  waves  into 
two  characteristic  magneto-ionic  modes  is  an  important  one  of  these. 
Depending  upon  the  geometry  of  the  propagation  path  and  the  geomagnetic 
field,  these  modes  can  differ  by  as  much  as  a MHz  in  the  MOF  to  a given 
range.  The  difference  in  radar  delay  associated  with  magneto-ionic  multi- 
moding is  of  order  a few  microseconds  for  typical  OTH  radar  propagation, 
so  it  remains  unresolved  unless  the  radar  has  a bandwidth  that  is  a sub- 
stantial fraction  of  a MHz.  If  nnresolvea,  this  multimoding  manifests 
itself  as  a fading  effect  observed  in  echoes  from  linearly  polarized 
scatterers . 

Other  forms  of  multimoding  have  been  invoked  to  explain  other 
forms  of  echo  fading  and  dispersion  in  time.  Usually, their  effect  is  only 
to  induce  signal  fading.  This  will  be  d.  cussed  further  in  a later 
section  where  the  use  of  fading  rate  as  an  identifying  characteristic  of 
propagation  mode  is  considered. 
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Figure  2-5.  Skywave  Backscatter  f-Scan  in  Which  Multiple-Hop- 

Propagated  F.ciioes  from  a Repeater  Appear  Prominently. 
The  upper  trace  is  caused  by  two-hop  propagation 
to  and  from  the  repeater.  The  middle  trace  is 
caused  by  tne-hop  propagation  in  one  direction 
and  two-hop  in  the  other. 
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2.2  Mode  Regions  and  Boundaries 

It  is  useful  in  the  ensuing  discussion  to  define  some  terms  and 
concepts  for  the  special  purpose  of  dealing  with  the  OTh  radar  mode  identifi- 
cation problem.  At  this  point  we  introduce  two;  these  are  the  mode's 
observable  regioi.  and  the  mode's  coverage  region.  Each  propagation  mode 
has  associated  with  it  one  of  each  kind  of  these  regions.  The  MOR  (Mode- 
Observable  Region) is  defined  in  the  radar  delay-frequency  plane  and  is 
the  locus  of  all  points  on  that  plane  at  which  a radar  echo  can  appear  as 
the  result  of  energy  that  has  propagated  via  that  particular  mode  and  has 
been  scattered  by  a target  at  the  earth's  surface.  The  latter  proviso 
can  and  must  be  made  in  order  to  limit  the  region  to  that  in  which  targets 
of  interest  can  occur.  Aircraft  fly  sufficiently  near  to  the  surface  to 
produce  echoes  only  in  the  node -observable  region,  while  meteors  and  space 
vehicles  can  produce  echoes  outside  the  region. 

The  second  term,  mode-coverage  region,  is  used  to  denote  a region 
in  the  frequency-distance  plane  which  represents  the  space  covered  by  the 
radar  at  any  frequency,  using  a particular  propagation  mode  in  its  operation. 
The  mode-coverage  region  can  be  derived  from  the  mode-observable  region  by 
applying  the  coordinate  conversion  process,  appropriate  to  that  mode,  to 
the  radar  delay  coordinate.  However,  because  the  relationship  of  radar 
delay  to  ground  range  differs  from  one  mode  to  another,  the  relationship 
among  the  various  mode-observable  regions  (for  different  propagation  modes) 
will  not  be  the  same  as  that  among  the  mode-coverage  regions.  For  example, 
an  overlapping  of  two  & ode-observable  regions  does  not  necessarily  imply 
a corresponding  overlapping  of  the  same  modes'  coverage  regions.  In  the 
discussion  thus  far,  the  radar  observable-angle  with  respect  to  the  array 
and  the  spatial  coordinate-bearing  from  the  radar  have  not  been  mentioned. 
However,  the  structure  of  the  ionosphere  and,  consequently,  the  required 
conversion  process  between  radar  delay  and  ground  range  can  vary  with 
azimuth.  Were  the  variation  rapid  enough,  it  would  be  necessary  to  define 
the  mode-observable  and  mode-coverage  regions  in  a three-dimensional 
space.  We  here  adopt  the  usual  practice  of  assuming  that  the  effect  of 
azimuthal  variation  can  be  handled  adequately  by  simply  redefining  the 
mode- observable  and  mode -coverage  regions  for  each  angle  from  the  antenna. 

2-11 


• • *-  i i-  • 


2.2.1  Types  of  Boundaries 

As  noted  above,  the  delay  time- frequency  plane  can  be  partitioned 
into  mode-observable  regions,  each  a characteristic  of  a particular  propa- 
gation mode.  The  overlapping  portions  are  defined  as  ambiguous  regions. 

The  location  of  the  regions' boundaries  in  the  delay  time -frequency  domain 
may  be  established  by  observation  or  by  derivation  from  a model  ionosphere. 

Segments  of  the  boundary  encompassing  the  regions  fall  into  two 
categories,  hard  or  soft.  Hard  boundaries  are  those  which  are  caused  by 
well-defined  geometrical  limiting  effects  such  as: 

a.  The  leading  edge  of  ground  backscatter. 

b.  The  caustic  focusing  region. 

c.  The  trace  of  the  minimum  propagation  elevation  angle 
in  the  delay  time- frequency  domain  for  the  lowest 
ionosphere  layer. 

d.  The  focusing  at  maximum  range  ( Croft.  1967). 

The  hard  boundaries  can  be  directly  computed  by  ray  tracing  techniques. 

With  the  exception  of  the  minimum  elevation  angle  of  the  lowest  layer,  the 
remaining  boundaries  are  often  readily  observed  on  backscatter  records  as 
enhancements  of  echo  strength.  Figure  2-6  illustrates  an  example  of  hard 
boundaries.  An  ambiguity  region  is  delineated  by  the  leading  edge  of 
the  E-mode  backscatter  and  the  focusing  at  the  maximum  range  of  the  F layer. 

The  soft  boundaries  are  those  boundaries  at  which  propagation  is 
geometrically  possible, but  at  which  the  noise  or  path  loss  has  degraded 
the  signal-to-noise  ratio  to  a degree  such  that  the  probability  of  detection 
is  voo  small  to  be  useful.  The  maximum  delay  time  trace  on  a backscatter 
record  corresponding  to  the  minimum  elevation  angle  for  an  FI  or  F2  layer 
(E  layer  present)  is  an  example  of  soft  boundary.  System  effects  such  as 
antenna  gain,  equipment  loss  or  gain  and  ionospheric  loss  mechanisms  modify 
the  trace  of  a soft  boundary.  Figure  2-7  illustrates  an  example  of  a soft 
boundary.  Figure  2-8  chows  hard  and  soft  boundaries  established  by  ray- 
tracing in  a complex  ionosphere  composed  of  three  superimposed  quasi- 
parabolic layers.  The  leading  edge  of  the  E,  F,  and  F2  layers  are  computed 
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Figure  2-6.  Backscatter  Trace  Illustrating  Hard  Boundaries 
and  an  Ambiguity  Region 
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Figure  2-8.  Ambiguous  Regions  Established  by  Hard  and  Soft  Mode  Boundaries 

Antennas  with  little  gain  near  the  horizon  will  produce  a soft  boundary  at, 
for  example,  the  longer  time  delay  boundary  of  the  E layer  mode-observable 
region.  This  effect  was  not  considered  in  making  this  diagram,  so  this 
boundary  is  indicated  to  be  a hard  one.  - 


by  finding  the  corresponding  delay  time  trace  where  the  condition 


where  T = delay  time 

a = elevation  angle 

holds.  The  zero  degree  ray  establishes  a segment  of  the  E mode  boundary. 
The  trailing  edge  soft  boundaries  may  be  established  with  the  aid  of 
system  sensitivity  contours  such  as  that  shown  in  the  figure.  The  constant 
loss  contours  are  superimposed  on  the  E and  F regions.  These  loss  contours 
follow  the  general  shape  of  the  leading  edge  of  the  clutter,  and  the  loss 
decreases  to  a minimum  value  near  the  leading  edge.  This  decrease  in  echo 
strength  corresponds  to  the  fadeout  at  lower  frequencies  which  occurs  in 
Figure  2-7.  The  soft  boundary  is  the  locus  of  points  in  the  frequency- 
radar  delay  plane  at  which  the  signal- to-noise  ratio  is  only  marginally 
adequate  for  useful  observation. 

In  the  above  example,  the  soft  boundary  at  the  low  frequency  (or 
long  delay)  edge  of  the  mode-observable  regions  was  imposed  by  propagation 
path  loss  that  reduced  signal  strength  to  a degree  at  which  it  became 
unusable.  We  may  now  consider  another  example,  in  which  a soft  boundary 
is  produced  by  another  mechanism. 

Figure  2-9  is  a computed  skywave  backscatter  trace  for  a sumner 
daytime  sounding  in  the  southward  direction  from  Whitehouse.  The  caustic- 
focus  traces  of  the  FI  and  F2  layers  are  nearly  coincident,  while  that 
of  the  E layer  lies  about  6 MHz  higher  in  frequency  than  do  the  F layer 
traces.  Traces  for  what  would  be  discrete  targets  at  range  of  2354  km  and 
1063  km  are  also  shown.  This  figure  indicates  the  high  frequency  "hard" 
boundaries  imposed  by  the  geometrical  optics  of  the  situation,  but  it 
does  not  indicate  other  boundaries  of  the  mode-observable  regions. 

Figure  2-10  shows  a plot  in  the  elevation  angle-radar  delay  (a-x)  plane 
for  rays  of  a 14-MHz  signal.  This  diagram  indicates  the  obscuring  effect 
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Figure  2-10.  Traces  in  the  Radar  Delay-Elevation  Angle  Plane  of  Rays 
Representing  14-MHz  Propagation  via  E,  and  F-  1-aycrs 
to  the  Earth's  Surface.  Refocusing  is  indicated  by  decreased 
slope,  focusing  by  increased  slope. 
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of  underlying  layers.  In  particular,  note  that  the  lower  portions  of 
the  curves  for  the  FI  and  F2  layers  enter  a cusp,  near  12  ms  and  15.5  ms, 
respectively,  that  represents  the  transition  to  the  underlying  layer. 

As  the  cusp  is  approached,  the  curve  approaches  the  horizontal.  This 
is  indicative  of  defocusing.  Thus,  the  strength  of  ground  clutter 
echoes  corresponding  to  these  points  on  the  curve  can  be  expected  to  be 
weak.  For  example,  for  the  F2  layer,  a strong  focusing  is  indicated 
at  about  13  ms,  and  as  the  range  increases  beyond  13  ms,  the  low  angle 
rays  become  increasingly  spread  and  would  fade  completely  by  the  time 
a 20  ms  range  is  reached.  This  is  a good  example  of  a soft  boundary 
attributable  to  the  obscuring  effect  of  an  underlying  layer.  The 
boundary  is  "soft"  because  defocusing  leads  to  a gradual  diminution 
of  echo  strength  with  increased  range.  The  dashed  line  in  Figure  2-9 
indicates  how  this  soft  boundary  would  limit  the  F2  mode  observable 
region.  The  precise  position  of  this  boundary  would  vary  with  those 
factors  influencing  the  radar's  sensitivity  such  as  path  loss,  noise 
level,  and  target  cross  section. 

In  some  cases,  the  maximum  delay  of  a mode's  observable  region 
occurs  at  a hard  boundary- -one  at  which  focusing  occurs.  The  maximum 
range  (or  delay)  focusing  has  been  studied  by  Crgfl  (1967)  who  has  shown 
that  its  occurrence  depends  upon  the  distribution  of  electrons  between 
two  ionospheric  layers.  Examples  of  this  can  frequently  be  observed  in 
the  frequency  scans  of  an  OTH  radar;  Figure  2-11  is  a good  example  of  this. 
The  time  axis  of  the  oblique  ionogram  is  too  compressed  to  indicate  the 
details  of  the  repeater-generated  oblique  ionograms  in  the  vicinity  of 
the  trace  of  the  maximum  range-focused  clutter.  Croft  has  treated  the 
matter  in  some  detail  and  has  found  that  an  oblique  ionogram  trace  will 
undergo  an  S-shaped  turn  in  making  the  transition  from  the  low  ray  trace 
of  an  overlying  layer  to,  ultimately,  the  low  ray  trace  of  the  underlying 
layer.  The  extension  of  the  repeater  trace  beyond  the  trace  of  the  maxi- 
mum range-focused  clutter  is  caused  by  low-ray  propagation  supported  by 
the  underlying  layer,  so  it  can  be  seen  to  occur  at  a slightly  lesser  delay 
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From  the  above  discussion,  we  may  summarize  the  following  points: 

(1)  The  leading  edge  (in  radar  time  delay)  of  a mode-observable 
region  is  usually  delineated  by  the  hard  boundary  produced 
by  the  familiar  focussing  effect. 

(2)  The  trailing  edge  of  a mode -observable  region  may  be  either 
a "hard"  or  a "soft"  boundary,  depending  upon  the  structure 
of  the  layer  and  upon  other  contingencies.  A soft  boundary’s 
position  will  vary  with  path  loss,  noise  level,  and  the 
radar  cross  section  of  the  particular  target.  Consequently, 
its  position  can  be  predicted  and  used  in  mode  identification 
decisions  only  with  careful  consideration  of  statistical 
detection  theory. 

2.2.2  Occurrence  and  Extent  of  Ambiguous  Regions 

In  this  section,  some  examples  of  the  extent  of  ambiguous  regions 
in  the  frequency- radar  delay  plane  are  given.  A comprehensive  analysis 
of  the  occurrence  of  ambiguities  is  not  intended;  only  a demonstration  that 
the  problem  occurs  a significant  fraction  of  the  time. 

Nighttime  conditions  at  temperate  latitudes  represent  the  simplest 
category  of  propagation  with  which  the  OTH  radarist  must  contend.  The 
underlying  FI  and  E layers  are  controlled  by  sunlight,  and  their  effect 
disappears  shortly  after  sunset.  The  F2  layer,  although  also  produced  by 
insolation,  responds  with  a significant  time  lag  (on  the  order  of  four  hours! • 
Consequently,  residual  ionization  prevails  through  much  of  the  night,  so  the 
ionosphere  has  the  form  of  a single  layer.  Figure  242  is  an  example  of  a 
sky  wave  backscatter  trace  obtained  under  these  simple  layer  conditions. 

Note  that  the  trace  of  repeater  returns  at  12  ms  extends  continuously 


from  the  minimum  frequency  used  by  the  radar,  to  the  boundary  of  the  mode- 
observable  region  delineated  by  the  minimum-range  focussed  clutter.  Were 
an  operating  frequency  near  13  MHz  used,  the  radar  would  cover  an  interval 
of  radar  delays  extending  from  about  8 ms  to  greater  than  17  ms.  The  latter 
figure  would  be  a soft  boundary  and  cannot  be  precisely  specified.  Figure  2-1 
shows  a similar  trace  with  improved  time  resolution.  Note  that  although 
the  pure  high -ray  trace  is  not  apparent,  there  appears  a trace  caused  by 
mixed  high -ray  and  low -ray  propagation.  The  mode- observable  region  for  this 
mixed-mode  extends  from  the  hard  boundary  at  the  trace  of  the  mi n imum -range  - 
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Ionosphere,  Typical  of  Winter  Nighttime  Conditions 


25  April  1971 
0229  Z 


Frequency  (Mils) 


Figure  2-13.  High  Resolution  Skywave  Backscatter  f-Scan  Made 
when  the  Ionosphere  Consisted  of  a Single  Layer. 

A region  of  ambiguity  occurs  because  of  high  ray- 
low  ray  mixed  propagation  within  3 MHz  of  the 
leading  edge  hard  boundary. 


focused  clutter  which  it  shares  with  the  mode  observable  region  (MOR)  of 
the  pure  low  ray  mode,  to  a "soft"  boundary  which  for  the  repeater  signal 
lies  approximately  3 MHz  below  the  frequency  of  the  hard  boundr  *v.  Measured 
on  the  radar  delay  axis,  the  soft  boundary  of  the  high- low  ray  MOR  lies  about 
2.5  ms  beyond  the  "hard"  minimum  delay  boundary.  This  high-low  ray  MOR  con- 
stitutes an  ambiguous  region  because  it  is  a subregion  of  the  low-low  ray  MOR. 
Consequently,  if  as  in  the  example  mentioned  earlier,  a 13-MHz  o;  ating  fre- 
que  .cy  is  used  for  the  radar,  the  interval  from  8 to  10.5  ms  would  be  ambiguous 
and  the  region  10.5  to  12  ms  or  more  would  not  be.  However,  the  trailing 
boundary  of  the  ambiguous  region  is  soft,  so  for  targets  with  cross  sections 
smaller  than  that  simulated  by  the  repeater,  the  ambiguous  region  would  be  less 
extensive.  Figure  2-14  shows  a similar  repeater-generated  oblique  sounding 
record  over  the  same  Whitehouse-to-Lincoln  path,  but  this  time  for  late  after- 
noon conditions.  Again,  the  high-low  ray  mode  extends  about  3 MHz  below  the 
hard  boundary  formed  by  focusing  at  the  mode’s  leading  edge.  The  extent  of 
the  ambiguous  region  is  about  3 ms  in  this  case.  There  is  some  indication  of 
a high-high  ray  mode  trace,  but  its  strength  is  significantly  less  than  that 
of  the  other  modes. 

A set  of  measurements  by  Tornatore  (1972)  indicates 
that  the  low-high  mode  would  be  weaker  than  the  low- low  mode  by  about  3 db, 
one  MHz  below  the  maximum  frequency  limit,  and  about  6 db  at  2 MHz  below  the 
hard  boundary  (MUF) . Thus,  a radar  sensitive  enough  to  have  a 3-db  margin 
of  detection  for  targets  of  interest  would  have  an  ambiguous  region  extend- 
ing about  1 MHz  below  the  maximum  frequency , or  a radar  delay  of  about  one  ms 
beyond  the  near-range  focused  clutter  trace.  This  kind  of  ambiguity  can 
almost  always  be  expected  even  with  the  simplest  ionospheric  structure. 

Also,  multimoding  and  ambiguity  can  occur  under  single  layer 
conditions  because  of  multiple-hop  modes.  Figure  2-15  is  an  example  of  such 
a skywave  backscatter  record.  A repeater -generated  trace  produced  by  single 
mode  propagation  occurs  at  12.2  ms  radar  delay  and  joins  with  the  high-low 
ray  trace  at  the  focused  clutter  trace.  A second  repeater  trace,  considerably 
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Figure  2-15.  Skywave  Backscatter  f-Scan  in  Which  an  Ambiguity 
Occurs  Because  of  Multiple-Hop  Propagation 


2-26 


-§&* 


■*:  t-*W^>*^5$ 


dispersed,  appears  at  ranges  from  12.9  to  13.7  ms.  Barely  perceptible  is 
yet  a third  trace  extending  in  radar  delay  from  13.5  to  beyond  14  me.  This 
latter  trace  results  from  propagation  to  and  from  the  repeater  via  a two-hop 
mode,  while  the  stronger  trace  is  the  result  of  a mixed  mode  (e.g.,  propa- 
gation to  the  target  via  a one-hop  mode  and  return  via  a two-hop  mode) . In 
this  case  the  observable  region  for  the  mixed  one-hop,  two-hop  mode  has  a 
high  frequency  boundary  that  is  about  9 MHz  below  that  of  the  one-hop  mode. 

This  difference  in  frequency  will  vary  with  the  range  of  the  target  (or 
repeater)  being  considered  and  with  the  details  of  the  layer's  electron 
density  distribution 

For  normal,  effective  heights  of  the  F layer  (under  400  km),  the 
maximum  frequency  at  which  two-hop  propagation  is  supported  to  any  given  range 
is  usually  less  than  75%  of  the  frequency  at  which  one-hop  propagation  is 
supported.  Under  nighttime  conditions,  when  the  F layer  critical  frequency 
can  be  as  low  as  3 MHz,  the  4000  km -range  MUF  can  be  as  small  as  10  to  12  MHz, 
resulting  in  a frequency  interval  between  the  MUF  of  the  one-hop  and  two- 
hop  modes  of  approximately  3 MHz.  The  difference  in  frequency  between  the 
leading  boundary  of  the  one-hop  M0R  and  the  mixed  ore-hop,  two-hop  M0R  would 
be  nearly  as  great,  perhaps  about  2.5  MHz.  It  is  conceivable  that  under 
these  conditions  a well  developed  one-hop  high-ray  trace  in  combination  with 
such  a minimal  frequency  difference  between  the  leading  boundaries  of  the  one- 
hop  and  two-hop  mixed  mode  MOR's  could  preclude  the  possibility  of  choosing  an 
operating  frequency  at  which  the  radar  delay  of  interest  would  be  free  of 
ambiguity;  however,  a brief  search  through  available  skywave  backscatter  records 
failed  to  discover  any  such  examples.  A number  of  examples  were  found  in  which 
a region  of  one-hop,  low-ray  F-mode  propagation  occurred  unambiguously  over 
a 3-MHz  frequency  interval . 

To  summarize  this  discussion  of  the  simple,  single -layer  ionosphere: 

(1)  Such  conditions  occur  frequently  at  night. 

(2)  The  one-hop,  low  ray  MOR  (M0R-1FL)  extends  from  the 
frequencies  of  7 MHz  or  less  to  the  hard  boundary 
delineated  by  the  focused  clutter  trace. 

(3)  MOR- IF!,  is  ambiguous  with  the  one-hop  high-ray  MOR 
(M0R-1FH)  in  a region  extending  1 to  2 MHz  below  the 
one-hop  F focussed  clutter  trace.  The  low  frequency 
boundary  of  the  ambiguous  region  is  soft. 
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(4)  The  M0R-1FL  is  ambiguous  with  the  mixed  one-hop, 
two-hop  MOR  (M0R-1F2F)  at  frequencies  less  than 
60%  to  75%  of  the  1FL  mode  HOF  for  any  range. 

(5)  The  high  frequency  boundary  of  the  M0R-1F2F  is  a 
hard  boundary,  but  it  is  usually  not  discernible 
in  the  skywave  backscatter  records  as  a trace  of 
focused  clutter. 

Figure  2-16  diagrams  the  MOR's  discussed  above.  At  least  a short 
interval  of  radar  delay  will  usually  be  free  of  ambiguities  when  such  con- 
ditions prevail. 

The  next,  more  complicated,  circumstance  is  that  in  which  the 
F layer  is  accompanied  by  an  underlying  layer.  This  may  be  either  the  E layer 
at  100  km  or  the  FI  layer,  which  occurs  near  200  km.  Of  course,  both  of  these 
underlying  layers  exist  simultaneously,  but  for  the  moment  we  will  consider 
the  effects  of  only  one. 

A sufficiently  dense  underlying  layer  could  entirely  obscure 
the  F layer,  so  no  trace  indicative  of  the  F layer  would  appear  on  the  record. 
This  does  not  happen  normally  but  may  happen  in  cases  of  intense  sporadic  E. 

A more  common  occurrence  is  for  the  underlying  layer  to  obscure  all  but  those 
backscatter  echoes  that  occur  in  the  vicinity  of  the  leading  focused  clutter 
trace.  An  example  of  this  is  shown  in  Figure  2-17.  An  expanded  trace  that 
better  shows  the  limiting  of  the  backscattered  echoes  to  within  a MHz  of  the 
leading  edge  is  provided.  As  noted  earlier,  the  region  of  the  one-hop  F MOR 
within  a MHz  of  the  leading  edge  is  subject  to  ambiguity  involving  the  high 
and  low  rays,  so  this  represents  a case  where  potential  ambiguity  is  unavoid- 
able. The  portion  of  the  E region  MOR  at  frequencies  higher  than  the  F 
layer's  focused  clutter  trace  is  unambiguous  and  could  be  used  by  the 
radarist.  However,  the  ranges  of  the  scatterers  observed  via  this  mode  are 
substantially  less  than  those  observed  via  the  F mode.  The  low  frequency 
boundary  of  the  E region  MOR  is  a soft  one  caused  by  absorption  of  the 
propagated  waves  in  the  I)  layer.  Figure  2-18  is  a schematic  diagram  showing  the 
distribution  of  the  mode-observable  regions  in  the  radar  delay-frequency  plane. 
The  high  frequency  edge  of  the  F layer  MOR  is  shown  as  a hard  boundary,  but 
the  oblique  ionogram  indicates  a raggedness  that  makes  the  position  of  the 
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Figure  2-17.  Skywave  Backscatter  f-Scan  Showing  How  Echoes  can  be  Restricted  to  the  Immediate 
Vicinity  of  the  Leading  Edge  of  an  M0R.  Figure  2-17(b)  is  a high-resolution 
presentation  of  a portion  of  Figure  2-1 7fa). 


boundary  uncer*  n.  This  may  be  the  effect  of  irregular  obscuration 
by  the  underlying  sporadic  F.  and  to  various  scattering  effects  caused  by 
the  same  layer. 

In  the  preceding  discussion,  the  underlying  layer  had  an  MOR 
that  extended  to  higher  frequencies  than  did  the  F layer  MOR.  Ke  next  con- 
sider two  cases  where  the  F layers  MOR  extends  to  frequencies  higher  than 
does  the  MOR  of  the  underlying  layer.  These  cases  were  discussed  in  the 
earlier  section  on  MOR  boundaries  and  to  much  greater  length  by  Croft  (1967). 
Here,  we  simply  note  that  one  or  the  other  of  these  cases  will  occur  during 
winter  daytime  conditions,  when  the  F layer  critical  frequency  is  anomalously 
high,  and  the  E layer  critical  frequency  is  somewhat  reduced  by  the  sun's 
greater  zenith  angle.  Figure  2-11  showed  an  example  in  which  the  low  frequency 
edge  of  MOR-1FL  occurs  at  a focused  clutter  trace,  while  the  example  in 
Figure  2-  Jy  is  of  a MOR-1FL  that  is  softly  bounded,  so  the  repeater  trace 
fades  away  at  its  low  frequency  limit. 

In  either  case  the  frequency  extent  of  the  MOR- 1 FI.  increases  with 
radar  delay,  being  about  4 Ml  I z at  a radar  delay  of  7 ms  and  increasing  to 
S MHz  at  12  ms  in  Figures  2-1  i and  2-19  . Also,  in  each  example,  the  low 
frequency  portion  of  the  MOR-IFL  intersects  the  MOR- IE  and  M0R-1E,1F,  indicating 
a region  of  mode  ambiguity.  The  ambiguous  region  can  be  barely  discerned  as 
a region  of  enhanced  clutter  in  the  radar  delay  interval  from  6 to  10  ms.  A 
more  clear  delineation  of  the  ambiguous  region  appears  in  Figure  2-20  where  there 
is  some  minimum  range  focusing  of  the  clutter  seen  via  the  IE  mode. 

Figure  2-21  is  a schematic  diagram  of  the  distribution  of  mode  ob- 
servable regions  under  these  conditions.  Note  th3t  there  is  considerable 
opportunity  for  ambiguity  in  the  vicinity  of  the  low  frequency  edge  of  the 
MOR- IF,  where  it  overlaps  the  MOR's  of  the  E layer  and  is  near  the  MOR  of 
the  two-hop  F mode,  as  well  as  the  MOR's  of  the  mixed  modes.  However,  it 
may  also  be  noted  that  there  is  an  interval  of  several  MHz  in  which  the  MOR 
of  the  one-hop  F layer  low  ray  is  unambiguous. 
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Figure  2-19.  Skywave  Backscatter  f-Scan  Illustrating  a Softly 
Bounded  Low  Frequency  Limit  to  the  Mode  Observabl 
Regions 
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Figure  2-20.  Skywave  Backscatter  f-Scan  Illustrating  an  Ambiguous 
Region  Where  M0R-1FL  Intersects  M0R-1E  and  M0R-1E,1F 
Minimum  range  focusing  of  the  IE  mode  is  apparent. 
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Figure  2-21.  Schematic  Diagram  of  Moiie  Observable  Regions  which  Ixi^t 

when  the  Higher  (l:)  Layer's  MOR  Extends  to  Higher  !•  requeue io 
than  does  the  Underlying  Layer's 
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The  three  classes  of  ionosphere  discussed  above  include  an 
important  fraction  of  the  situations  with  which  an  OTH  radar  must  deal, 
respectively  representing  nighttime  conditions,  summer  daytime  and  winter 
daytime.  An  experienced  HF  radarist  will  recognize  that  the  simple 
categorization  made  here  does  not  always  hold.  For  example,  during  summer 
nights  the  sporadic  E layer  may  introduce  additional  modes  in  the  manner  of 
the  daytime  E layer,  and  during  ionospheric  storms  the  F layer  critical 
frequencies  may  be  reduced  by  15%  to  several  tens  of  percent  (Matsushita, 

1959)  depending  upon  the  geomagnetic  latitude  of  the  reflection  point. 

Also,  there  are  the  transitions  between  daytime  and  nighttime  conditions. 
During  what  we  have  here  proferred  as  typical  winter  conditions,  no  serious 
problem  arises  because  the  M0R-1FL  for  both  daytime  and  nighttime  conditions 
extends  to  higher  frequencies  than  do  those  of  the  underlying  layers.  How- 
ever, during  summer  conditions,  a transition  must  take  place  wherein  the 
high  frequency  boundary  of  the  E layer  MOR  must  move  across  that  of  the  F 
layer.  At  ranges  beyond  those  at  which  one-hop  E can  occur,  there  is  still  no 
problem.  The  only  effect  is  the  one  of  obscuration,  which  does  not  occur 
until  the  E layer  MOR  at  any  given  radar  delay  extends  to  appreciably  higher 
frequencies  than  does  the  F layer's.  However,  3t  E layer  ranges  (out  to 
about  1200  km  or  about  ?0  nu  radar  delay)  a well  developed  ambiguity  occurs. 
Figure  2-22  diagrams  the  MOR's  that  occur  during  such  a transition.  During  th 
hour  following  1800,  the  E layer  MUF  drops  rapidly  and  by  1900  lies 
a few  MHz  lower  in  frequency  than  the  F-layer  MUF. 

This  completes  the  discussion  of  the  occurrence  and  extent  of 
mode  ambiguity.  We  have  noted  that,  provided  one  is  free  to  select  the  radar'; 
operating  frequency  solely  on  the  criterion  of  nonambiguity  of  the  radar 
echoes  received,  there  usually  can  be  found,  for  any  radar  delay,  a few  MHz 
interval  free  from  ambiguities.  An  exception  to  this  sometimes  occurs  during 
summer  daytime  conditions  when  obscuration  by  underlying  layers  may  limit  the 
F layer  MOR  to  a band  extending  only  a MHz  or  so  below  the  focussed  minimum- 
range  clutter  trace.  Under  such  conditions,  propagation  utilizing  the  high- 
ray  mode  on  the  outward  or  inward  path  produces  a delayed  trace  that  has  a 
sifnal  strength  on  the  order  of  7>  db  '.''ss  than  that  of  the  echo  propagated 
entirely  by  the  low  ray  mode. 
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Mode  ambiguity  occurs  also  at  short  ranges  at  which  E layer 
propagation  can  occur.  A particular  problem  is  posed  by  periods  of  day/ 
night  transition,  when  the  leading  edge  traces  cross. 
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TESTS  FOR  MODE  IDENTIFICATION  INVOLVING  BOUNDARIES 


In  this  section  it  is  assumed  that  individual  detections  by  an 
OTH  radar  are  first  processed  and  assembled  into  tracks.  In  the  present 
context,  a track  is  a subset  of  the  detections,  which,  on  the  basis  of 
proximity  to  one  another  and  consistency  of  doppler  shift,  range  rate, 
and  direction,  appear  to  be  the  results  of  multiple  observations  of  a 
single  target.  If  the  tracking  is  done  in  the  space  of  radar  observ- 
ables (i.e.,  frequency,  radar  delay,  and  antenna  angle),  a single 
target  that  is  illuminated  by  multiple  modes  may  give  rise  to  multiple 
tracks.  If,  on  the  other  hand,  the  tracking  is  to  be  done  in  the  form 
of  geographic (or  other  true  space)  coordinates,  it  is  necessary  to  per- 
form coordinate  conversion  on  a detection-by-detection  basis.  In  this 
case,  if  a detection  happens  to  lie  within  an  ambiguous  region  of  the  T-f 
plane,  im.  xtiple  hypotheses  would  have  to  be  carried  regarding  the  operative 
mode  of  propagation,  resulting  in  multiple  hypothetical  positions  for  the 
detected  target,  and  consequently,  multiple  tracks.  Thus,  in  either 
event  the  effect  of  multiple  propagation  modes  is  to  raise  the  possibility 
of  there  being  more  radar  tracks  than  there  are  targets.  For  the  present 
discussion,  we  assume  that  the  tracking  is  to  be  done  in  the  x-f  plane. 

Also  we  begin  with  the  discussion  of  the  more  general  case  of  a frequency- 
hopping radar.  The  problem  of  mode  identification  for  a fixed-frequency 
radar  is  covered  In  later  portions  of  this  report. 


3. 1 Tracks  in  the  f-x  Plane 

If  a fixed  target  is  observed  by  a frequency-hopping  radar,  the 
individual  observations  will  lie  on  a trace  that  is  equivalent  to  an 
obliquely  sounded  ionogram  made  over  the  path.  The  detections  of  an 
approaching  (or  receding)  target  would  form  a cloud  of  points  that  would 
be  bounded  by  the  oblique  ionogram  traces  for  paths  to  the  near  and  far 
(with  respect  to  the  radar)  points  of  the  target  during  the  period  of 
observation.  (See  Figure  3-1.) 


figure  3-1.  Occurrence  of  Detections  in  the  x-f  Plane  When  an 
Approaching  Target  is  Observed  with  a l-’requency- 
liopping  Radar 
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If  the  radar's  range  resolution  (bandwidth)  is  not  great  and  if 
the  period  of  observation  is  sufficiently  small,  the  displacement  during 
the  course  of  the  observations  may  be  small  compared  to  the  radar's 
resolution,  and  the  observations  may  all  fall  on  a single  oblique  ionogram 
trace  (although  given  sufficiently  great  signal -to-noise  ratio,  the  dis- 
placements could  be  measured).  Assuming  10-second  dwells  and  a 900  km/hr. 
aircraft,  the  displacement  would  be  about  2.5  km  during  a dwell  and 
would  be  discernible  from  one  observation  to  the  next,  only  if  the  radar's 
bandwidth  exceeded  60  kHz.  With  a 20-kHz  radar,  an  observation  at  any 
time  would  have  a radar  delay  discemibly  different  from  those  made  more 
than  two  dwells  earlier  or  later,  while  a radar  with  a 250-kHz  bandwidth 
would  see  an  appreciable  change  in  range  from  one  dwell  to  the  next. 

Also,  a wide  bandwidth  system  (-250  kHz)  would  experience  the  effects  of 
dispersion  in  the  ionosphere  so  that  observations  of  a target  made  at 
frequencies  near  the  MUF  would  be  at  discernibly  greater  radar  delays 
than  would  those  made  well  below  the  MUF. 


The  above  properties  are  of  importance  to  the  tracker,  which  must 
take  account  of  them  either  by  correct  prediction  of,  or  by  being  tolerant 
of,  the  effects. 


3.2  Tests  for  Mode  Identity 

The  mode  identification  tests  to  be  discussed  in  this  section  all 
involve  the  behavior  of  the  signatures  at  the  boundaries  of  the  mode 
observable  regions.  The  basis  for  these  tests  is  that  a track  must  be 
entirely  confined  within  the  mode  observable  region  (MOR)  for  the  opera- 
tive mode.  If  a track  includes  points  outside  of  a particular  MOR,  then 
one  may  conclude  that  the  corresponding  mode  was  not  the  mode  that  pro- 
vided the  detections  which  made  up  the  track.  Conversely,  if  all  of 
the  detections  occurred  within  a particular  MOR,  then  one  may  be  justi- 
fied in  identifying  the  propagation  mode  of  the  detections  as  that 
corresponding  to  the  MOR.  However,  now  there  must  be  a certain  amount 
of  reservation  in  making  the  identification.  For  example,  the  set  of 
detections  may  also  all  be  within  a second  MOR  as  well;  i.e.,  the 


detections  may  all  lie  within  an  ambiguous  region.  In  order  to  resolve 
such  an  ambiguity,  observations  must  be  made  such  that  a detection  occurs 
(or  is  observed  not  to  occur)  within  an  unambiguous  region,  or  within 
some  combination  of  ambiguous  regions  which  provides  an  unambiguous 
result. 

Another  concept  is  that  the  credibility  of  an  unambiguous 
detection,  or  of  an  unambiguously  missed  detection,  must  be  considered. 

An  OTH  radar  often  operates  in  a noisy  environment  and  is  subject  to 
false  alarming.  Also,  the  signals  fade  because  of  various  propagation 
effects,  so  missed  detections  are  to  be  expected.  The  relative  prob- 
ability of  these  two  kinds  of  error  are  determined  by  the  thresholds 
employed  in  making  detection  decisions.  Clearly,  a few  detections 
outside  of  an  MOR  should  not  rule  out  the  associated  mode  if  the  likeli- 
hood of  such  detections  occurring  erroneously  is  large.  Similarly,  the 
absence  of  detections  outside  of  an  MOR  cannot  be  regarded  as  a con- 
clusive indication  of  mode  if  there  is  an  appreciable  probability  that 
the  detections  could  have  been  missed  by  chance.  Thus,  mode  identifica- 
tion by  boundary  crossing  involves  considerations  of  prolabilities,  even 
in  the  idealized  case  where  the  boundaries  are  well  defined  and  pre- 
dictable. At  soft  boundaries,  where  the  probability  of  detection  dimin- 
ishes gradually  with  frequency  or  radar  delay,  the  need  for  probabilistic 
arguments  is  all  the  more  acute. 

Having  discussed  some  of  the  general  principles  involved  in  this 
method  of  mode  identification,  we  now  consider  some  particular  examples. 
The  first  example  will  involve  the  overlapping  of  two  mode  observable 
regions  to  form  an  ambiguous  region.  Figure  3-2  indicates  how  the  MOR's 
might  be  related.  Two  traces  of  target  detections  are 
shown.  The  upper  trace  extends  beyond  the  high  frequency  boundary 
of  the  MOR  of  Mode  B.  Because  these  detections  cannot  have  occurred 
via  Mode  B,  it  may  be  concluded  that  Mode  A was  the  operative  mode.  On 
the  other  hand,  the  lower  trace  terminates  at  the  boundary.  Had  these 
detections  occurred  via  Mode  A,  then  they  would  not  be  expected  to  be 
influenced  by  the  Mode  B boundary.  Therefore,  we  may  conclude  that 
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Figure  3-2.  Region  of  Ambiguous  Detections  Procured 

by  the  Overlapping  of  Two  .lode  Observable 
Regions. 
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these  detections  were  made  via  Mode  B.  However,  we  again  emphasize  that 
the  occurrence  and  non-occurrence  of  the  detections  is,  to  some  degree, 
fortuitous,  and  probabilistic  considerations  must  be  taken  to  ascertain 
a trace's  extension  into  the  unambiguous  portion  of  the  A mode's  MOR. 

As  the  next  example,  consider  an  ambiguous  region  produced  by  the 
overlapping  of  three  MOR's.  In  principle,  this  could  occur  whenever  two 
propagation  modes  illuminate  a target  or  region  of  the  earth's  surface 
because  a mixed  (out  via  A,  back  via  B)  propagation  mode  would  occur  along 
with  the  two  unmixed  modes.  Figure  3-3  illustrates  the  MOR's  overlapping 
in  this  way.  Seven  traces  are  shown.  Traces  1,  2,  and  3 each  cross  two 
boundaries  and  can  all  be  identified  readily.  Trace  1 and  Trace  3 extend 
into  unambiguous  portions  of  the  mode-C  MOR,  and  so  they  can  be  identified 
as  being  the  result  of  propagation  via  mode  C.  Trace  2 extends  on  the  left 
into  a region  in  which  B-mode  detections  cannot  occur  and  on  the  right 
into  a region  in  which  A-mode  detections  cannot  occur.  Therefore, 
we  may  conclude  that  the  detection  was  effected  via  the  C-mode.  In 
somewhat  similar  fashion,  truncation  of  the  traces  at  a boundary  world  have 
permitted  other,  equally  definitive  conclusions  to  be  drawn.  Thus,  had 
Trace  1 terminated  at  the  left  (low  frequency)  boundary  of  the  A MOR,  so 
that  it  did  not  extend  into  the  unambiguous  region  of  the  C MOR,  we  would 
have  concluded  that  A was  the  operative  mode.  Similar  reasoning  would 
have  applied  to  Traces  2 and  3 had  they  been  observed  to  terminate  at 
an  MOR  boundary  rather  than  extending  beyond  it. 

Traces  4,  5,  6,  and  7 are  examples  of  traces  which  cross  a single 
boundary.  Traces  4 and  7 extend  into  unambiguous  subregions  of  the  C 
MOR  and  so  ..aist  have  been  effected  by  the  C mode.  However,  the  ambiguity 
of  Traces  5 and  6 cannot  be  resolved.  Each  of  these  traces  extends 
from  a region  of  three-fold  ambiguity  into  a region  of  two- fold  ambiguity, 
so  the  two-fold  ambiguity  remains.  Thus,  for  example,  the  extension  of 
Trace  6 beyond  the  ABC  ambiguous  region  into  the  BC  ambiguous  region  leaves 
the  choice  between  B-mode  and  C-mode  propagation  unresolved. 
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The  truncation  of  a trace  at  a boundary  yields  more  definitive 
information  because  it  uniquely  associates  the  trace  with  the  mode 
whose  MOR  is  bounded.  As  an  example,  if  Trace  5 were  limited  to  the  region 
ABC,  instead  of  extending  across  the  boundary  into  the  region  AC,  we 
could  conclude  that  the  trace  is  the  result  of  propagation  via  Mode  B 
at  whose  MOR  boundary  the  trace  terminated. 

In  summary,  the  preceding  discussion  has  identified  three  rules 
that  can  be  applied  to  permit  the  indentification  of  propagation  modes 
of  radar  detections  that  have  been  developed  in  the  radar-observables 
space  (T-f  plane).  They  are: 

(1)  Extension  of  the  trace  into  (or  total  inclusion  within) 

an  unambiguous  region  indicates  that  the  track  resulted  from 
propagation  via  the  unambiguous  regions  mode. 

(2)  Termination  of  a trace  at  the  boundary  of  an  MOR  indicates 
that  the  propagation  occurred  via  the  mode  of  the  MOR  in 
which  the  trace  is  contained.  In  applying  this  test,  it  is 
necessary  to  have  made  observations  adequate  to  have  dis- 
closed an  extension  of  the  trace  beyond  the  MOR's  boundary, 
had  it  existed. 

(3)  Extension  of  a trace  lying  within  a triply  ambiguous 
region  into  two  doubly  ambiguous  regions  permits  the  identi- 
fication of  the  propagation  mode  to  be  identified  as  that 
common  to  the  two  doubly  ambiguous  regions.  (It  is  assumed 
that  one  and  only  one  mode  recurs  in  both  of  the  doubly 
ambiguous  regions.) 

3 . 3 Probability  of  Detection  and  False  Detection 

The  boundary  tests  for  mode  identification  require  that  it  he 
determined  whether  or  not  a target-generated  trace  extends  beyond  the 
boundary  of  the  MOR  in  which  it  occurs.  Because  OTIi  radars  are  fre- 
quently operated  at  relatively  low  signal-to-noisc  ratios,  there  exists 
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some  chance  that  any  observation  will  give  an  erroneous  result, 
indicating  falsely  the  presence  of  a target  or  failing  to  indicate  a 
target  when  an  echo  is  indeed  present.  If  the  probability  of  such 
errors  is  greater  than  can  be  tolerated  for  the  probability  of  mis- 
identifying  the  mode  responsible  for  a trace,  it  is  necessary  to  defer  the 
mode  identification  decision  until  a sufficient  number  of  observations  can 
be  made  to  permit  a more  reliable  decision.  We  here  treat  the  probabili- 
ties of  these  decisions. 

One  approach  to  a tracking  algorithm  involves  a predictive  filter 
which,  on  the  the  basis  of  prior  observations,  extrapolates  the  observed 
target  positions  to  a position  at  which  it  is  anticipated  that  it  will 
be  detected  in  subsequent  observations.  As  additional  observations  are 
made,  the  new  information  is  used  to  improve  the  predictions  for  later 
observations.  The  logic  for  making  the  detections  in  such  a tracker  can 
take  various  forms.  For  example,  a region  of  acceptance  about  the 
predicted  point  can  be  defined,  and  within  that  region  the  cell  contain- 
ing the  largest  energy  can  be  regarded  as  the  new  detection.  Because 
one  cell  will  always  have  more  energy  than  the  others,  this  approach 
has  the  characteristic  of  making  a detection  each  time  and  would  be 
incapable  of  prouucing  a "no  detection"  result.  This  property  is  undesir- 
able for  our  present  purpose  because,  as  discussed  above,  the  radar's 
failure  to  detect  a tracked  target  in  any  given  observation  provides  a 
certain  kind  of  information  required  by  the  mode  identification  process. 

For  this  reason  it  is  important  that  the  detection  process  employ 
some  form  of  thresholding  process,  ('.ells  with  less  energy  than  the 
threshold  level  are  regarded  as  devoid  of  echoes.  If  two  or  more  cells 
exceed  the  threshold,  then  the  cell  with  the  largest  energy  can  be 
selected  as  the  detection.  A more  elaborate  scheme  permits  the  threshold 
to  be  increased  with  increased  distance  of  the  cell  from  the  point  at 
which  the  echo's  position  has  been  predicted.  In  this  way,  detections 
arc  favored  that  most  close!)  conform  to  the  prior  indications  of  the 
track. 


The  consequence  of  this  kind  of  thresholding  scheme  is  that  there 
will  be  some  probability  of  detection  when  echoes  are  present  and  some 
probability  of  false  alarm  when  they  are  absent.  There  are  many  con- 
tingencies to  the  problem,  aside  from  the  commonly  considered  one  of 
additive  system  noise  with,  say,  a Gaussian  distribution.  For  example, 
the  predictive  tracker  may  err  because  of  past  erroneous  detection;  or 
the  target  may  accelerate,  causing  the  echo  to  fail  outside  of  the 
region  of  acceptance,  or  so  far  from  the  predicted  position  that  its 
energy  fails  to  exceed  an  increased  threshold  level.  Also,  polariza- 
tion fading  effects  can  cause  the  echo  to  be  undetectable.  For  our 
present  purpose,  we  assume  that,  in  spite  of  these  several  contingencies, 
the  detection  process  may  be  regarded  simply  as  a Bernoulli  process, 
i.e.,  a random  process  with  two  possible  outcomes- -detection  or  failure 
to  detect--with  probabilities  1’j,  and  1 - PQ,  respectively.  In  the 
event  that  no  echo  is  present  in  the  cell,  we  may  replace  the  probability 
of  detection  with  the  probability  of  false  alarm  so  that  the  two  outcomes 
of  the  Bernoulli  process  become  Pp^  and  1 - Pp^.  Then  an  M of  N de- 
tection logic  can  be  applied  to  determine  the  presence  or  absence  of  a 
track  in  any  given  subset  of  the  data.  In  this  way,  a decision  can  be 
based  on  multiple  observations,  and  the  probabiltiy  of  error  can  he  re- 
duced correspondingly.  However,  the  certainty  with  which  decisions 
critical  to  mode  idenfieat ion  can  be  made  may  be  limited  because  they 
are  based  on  a subset  of  all  the  observations  of  a given  target--e.g  , 
those  in  which  detection  in  a region  of  reduced  mode  ambiguity  is  at- 
tempted. If  the  frequencies  of  the  radar  are  selected  so  that  there  an 
too  few  attempts  to  observe  echoes  from  the  target  in  a critical  region, 
a reliable  decision  is  precluded. 

It  is  desirable  that  the  mode  identification  be  made  with  a prob- 
ability of  correctness  comparable  to  that  required  of  the  decision  that 
a target  is  present  and  should  be  entered  onto  the  display  for  the  user 
of  the  radar.  However,  because  the  number  of  observations  available  in 
the  region  critical  to  mode  identification  is  always  a subset  of  the 
detections  used  for  the  decision  to  indicate  the  target  on  the  display. 


the  probability  of  correct  Jiv  i.Mon  is  always  less.  Because  mode 
identification  is  necessar>  before  the  radar’s  coordinates  can  be  con- 
verted to  geographic  coordinates,  one  must  compromise  by  either 

(1)  tolerating  a higher  probability  for  error  in 
mode  identification  than  for  target  indication 

(2)  deferring  display  of  the  target  until  a reliable 
mode  identification  can  be  made 

There  are  variations  of  these  possibilities  that  can  be  used  as 
well.  For  example,  targets  can  be  displayed  at  each  possible  one  of  the 
ranges  when  ambiguities  exist  or  van  be  displayed  prior  to  the  time  hhen 
a sufficiently  certain  mode  idem  1 1 k at  ion  can  be  made  if  the  target 
appears  to  be  within  a geographical  region  where  it  poses  a threat. 

Let  Pj  be  the  probabi l i ty  with  which  it  is  required  that  the 
mode  identification  he  correct  and  assume  that  the  decision  can  be  based 
on  the  detection,  or  failure  to  detect,  a target  in  a subset  of  N obserxa 
tions.  Also,  let  the  probability  of  detection  be  P()  and  the  probability 
of  false  alarm  be  P(  Then,  in  the  event  that  echoes  are  present  in  the 
critical  region,  the  probability  of  realizing  k detections  in  n attempts 
to  detect  is  b(k;n,P.,i  .in.  to 
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and  the  probability  thit  there  nil  I be  k^  or  more  detections  is  then 
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These  are  commonl)  used  .1  ist  • ii  .it  ions,  and  *an  be  found  tabulated.* 
Similarly,  the  probability  «.•!  experiencing  kQ  or  more  false  alarms  m n 
tries  is 
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Because  an  erroneous  decision  in  the  mode  identification  process 
has  a similar  effect,  whether  it  be  a false  alarm  or  a missed  detection,  the 
"cost"  of  the  two  kinds  of  errors  are  equal--the  mode  operative  in  the 
detection  is  misidentif ini,  and  the  target  is  interpreted  as  being  in  the 
wrong  place;  consequently,  wo  wish  to  select  values  of  kQ,  P^  and  PpA  such 
that  the  two  kinds  of  errors  have  equal  and  suitably  low  probabilities. 

Assuming  that  the  test.-,  of  regions  containing  echoes  occur  with  the  same 
frequency  as  tests  mud  • in  regions  where  the  echo  is  present,  the  probability 
of  correctly  identifying  the  mode  P^  may  be  written; 

P(k>ko;n,PD)  = P, 
and  P(k>ko;n,PpA)=  1 i'j 

In  general,  and  l*^  . parameters  determined  by  the  characteristics 
of  the  radar,  and  by  the  req  . rements  placed  on  its  probability  of  detection 
and  false  alarm.  Particularly  influencing  these  numbers  would  be  the  SNR 
of  the  radar  echoes,  the  number  of  cells  searched  by  the  tracker,  and  the 
relative  "cost"  of  false  alarms  and  missed  detections.  Note  that  the 
probability  of  false  alarm  employed  here,  Pp^,  applies  to  the  overall 
process  of  attempting  to  extend  the  track  with  the  current  observation. 

Thus,  if  the  tracking  ulgoiithm  attempts  to  extend  the  track  by  seeking  a 
detection  (with  uniform  thieshold)  in  any  one  of  n cells,  this  probability  is 
n Pj.^  where  PpA  in  t h • • mon-  lommonly  encountered  probability  of  a false 
detection  resulting  flow  a single  trial.  Typically,  an  OTH  radar's  tracker 
may  search  from  10  ru  loo  <■<.•!  Is  in  its  attempt  to  extend  the  tiack,  so  p, >A 
must  be  between  one  .uni  tw  - powers  of  10  less  than  P^. 

As  an  example,  w.ipi  o-i  that  the  signal-to-noi se  ratio  at  the  detector 
is  10  db.  Then,  it  ■dn.uid  be  possible  to  detect  targets  with  a probability  of 
0.0  while  experiencing  u one- trial  false  alarm  probability  of  about  It)'*'  If 
the  tracker  examini  i>"i  ^lls  m its  at  temp  to  extend  the  track,  tin  probabil- 
ity P...  of  erroneousl\  extending  a track  when  tho  signal  is  absent  can  be 
taken  as  10  '.  In  so.!,  u vu.>-,  if  as  few  as  2 detections  in  10  trial*  un- 
regarded as  indicative  of  a target,  the  probability  of  successfully  detecting 
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a trace  extension  is  P {S  ^ 2 |l0,  .6}  = 0.998  while  the  probability 
of  falsely  detecting  a path  extension  P { S 2 | 10,  .01}  , has  a value  less 
than  .0043,  which  means  the  probability  of  correctly  concluding  that  no 
echo  is  present  exceeds  .995.'. 

The  above  example  demonstrates  the  power  of  a compound  test  of 
the  m-out-of-n-trials  kind.  As  expected,  the  test  is  less  powerful  if 
fewer  trials  are  available  for  making  the  decision.  For  exan^le,  with 
the  same  probabilities  of  detection  and  false  alarm,  but  with  only  five 
trials  to  make  the  decision,  a criterion  of  1 detection  out  of  5 trials 
would  lead  to  a probability  of  correctly  extending  a trace  of  0.913, 
and  a probability  of  0.951  that  a correct  decision  not  to  extend  the 
trace  into  an  MOR  would  be  made. 

In  the  above  examples,  lere  only  the  incrementing  of  the  number 
m of  successes  of  a given  n trials  was  at  our  disposal,  the  resulting 
probabilities  of  the  two  kinds  of  errors  that  can  occur  differed  from 
one  another  somewhat.  Nonetheless,  the  best  choice  that  could  be  made 
in  each  of  the  cases  resulted  in  error  probabilities  that  would  be 
quite  satisfactory  for  many  purposes.  Figure  3-4  shows  various  combinations 
of  PpA  and  Pp  which  can  be  used  with  a corresponding  m out  of  n decision 
criterion  to  achieve  a probability  of  0.95  that  the  existence  of  a 
trace  in  a critical  MOR  is  correctly  determined.  Similar  curves  are 
readily  developed  for  other  values  of  Pj,  using  tables  of  binomial 
distributions. 

3.4  Considerations  for  Testing  at  Soft  Boundaries 

At  a "soft"  boundary  of  an  MOR,  the  signal-to-noise  ratio  gradually 
degrades  until  the  signal  becomes  effectively  non-detectable.  The  important 
difference  is  that  there  is  no  sharp  boundary  dividing  regions  where 
echoes  arriving  via  a particular  mode  can  be  expected  from  those  where 
they  are  unexpected.  However,  the  methods  described  above  can  be  applied 
if  critical  regions  are  correctly  defined.  Let  us  suppose  a two-out-of- 
five  decision  criterion  is  being  used  to  determine  trace  extensions  into 
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Figure  3-4.  PR  and  P[;^  Acquired  to  Yield  a Probability  Pj=0.95 

That  the  F.xistence  of  a Trace  Fxtension  into  a 
Critical  MOR  is  Correctly  Determined  for  Various 
fii  out  of  n Criteria. 
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critical  regions  and  that  the  desired  probability  of  correctly  making 
these  determinations  is  to  be  0.95.  Referring  to  Figure  3-4,  we  find 
that  the  probabilities  Pp  =0.66  and  PpA  = 0.07  are  required  to  achieve 
the  required  probability  of  a correct  decision.  A subset  of  five  obser- 
vations must  be  chosen  such  that  we  are  sure  of  attaining  this  degree 
of  certainty  regarding  the  outcome  of  the  individual  observation.  Let 
the  softly  bounded  mode  be  called  A,  and  let  B be  the  mode  with  which 
the  A mode  is  ambiguous.  The  critical  region  for  the  mode  identification 
test  must  be  chosen  so  that  if  the  echoes  have  arrived  via  mode  B,  they 
will  be  detected  by  each  observation  with  a probability  Pp  = 0.66.  The 
signal-to-noise  ratio  of  importance  here  is  that  of  the  B node,  which  for 
now  we  assume  to  be  uniform  over  the  B- mode- observable  region.  The 
combination  of  B mode  SNR  and  required  detection  probability  determines 
an  upper  limit  for  the  detection  threshold  Tp.  In  a simple  detection 
problem  a lower  limit  to  the  threshold  value  is  fixed  by  the  maximum 
false  alarm  rate  which  can  be  tolerated  when  no  signal  is  present.  A 
similar  consideration  applies  here,  except  that  the  false  alarm  event  is 
no  longer  caused  solely  by  noise,  but  rather  it  can  also  be  caused  by 
a signal  which  differs  from  what  would  be  a B-mode  signal,  only  in 
having  been  attenuated  by  the  mechanism  that  operates  to  bound  the  A- 
mode -observable  region.  In  order  to  avoid  a false  mode  identification 
using  two  of  five  decisions,  we  find  by  referring  to  Figure  3-4,  that  the 
probability  of  "false  alarming"  on  a single  observation  must  be  .07  or 
less. 

This  single  observation  involves  the  testing  of  all  points  within  some 
acceptance  region  determined  by  the  characteristics  of  the  tracker.  Let  us 
suppose  that  this  acceptance  region  comprises  N cells.  (N  will  usually  be 
somewhere  between  10  and  100.)  In  N-l  of  these  cells,  only  noise  is  present, 
while  in  one  of  the  cells  there  is  the  attenuated  signal  that  has  propa- 
gated via  the  A mode.  The  probability  of  making  a detection  based  on 
noise  only,  we  write  as  PpA  ^ ~ (N-l)  p^a  where  Pfa  is  the  familiar 
probability  of  false  detection  for  a single  threshold  comparison  of 
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the  signal  sample  from  a single  resolvable  cell  of  the  radar.  The 
probability  of  the  attenuated  echo  being  detected  in  the  remaining 
cell,  denoted  here  P'fa  is  just  the  probability  of  detection  of  a signal 
with  the  reduced  SNR.  Thus,  the  choice  of  the  threshold  for  this  decision 
is  determined  by  the  equation 

PFA  = CN-1)  pfa  ♦ pya 

where,  for  the  example  under  consideration,  PpA  = 0.07.  The  relative 
magnitudes  of  the  two  terms  will  vary  with  the  SNR  of  the  attenuated 
A-mode  signal.  Figure  3-5  indicates  how  the  reqqired  decision  threshold 
must  be  varied  with  the  anticipated  value  of  the  SNR  for  a given 
probability,  p'fa.  Figure  3-6  indicates  the  minimum  threshold  levels 
that  can  be  used  to  reject  echoes  of  the  specified  SNR  for  the  example. 

This  chart  was  derived  by  (1)  selecting  a threshold  level,  (2)  determining 
the  probability  of  false  alarm,  p^a>  in  24  trials  (assumed  tracker  search 
area),  (3)  subtracting  p^a  from  the  desired  overall  probability  of  "false'* 
detection,  PpA,  to  yield  residual  probability  p'^,  and  (4)  determining 
the  SNR  which  would  just  lead  to  the  residual  probability  p'^  at  the 
threshold  level  selected  in  (1) . 

The  relationships  among  threshold  levels,  probability  of  detection, 
and  SNR  were  determined  by  referring  to  Figure  3-5  and  to  graphs  of  the 
incomplete  Toronto  function  published  by  Marcum  (1960) . The  incomplete 
Toronto  function  T^  (1,  0,  /q)  is  equal  to  the  probability  of  not  detecting 

a single  pulse  signal  of  SNR  = q using  a threshold  v (measured  in  units  of  twice 
the  variance  of  the  noise).  Note  that  for  our  rather  typical  example,  the 
threshold  level  required  to  reject  a signal  of  SNR  = 3 (4.8  db)  with  the  de- 
sired probability  is  4 standard  deviations  above  the  mean.  At  this  threshold 
level  (v  = 8)  in  order  to  detect  unattenuated  (B-mode)  targets  with  the 
desired  probability  of  0.66,  it  is  necessary  for  their  SNR  to  be  9.5  (9.7  db) . 

In  the  absence  of  echoes,  the  probability  of  false  alarming  with  a 4a 
threshold  would  be  0.0003  on  a single  trial,  or  about  0.0075  in  a 
25-celi  search  area  of  the  tracker.  A threshold  level  of  3.5  standard 
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q (SNP  LINEAR) 

rc'hold  v Required  to  Rejc  ignals  with  SNR  of  q Derived  from  a Graph  of 
o incomplete  'ioronte  Runet  I^-v-  (1,  0,  ffi)  (Marcum,  1960) 


(STANDARD  DEVIATIONS  OF  NOISE) 


deviations  could  be  used  if  the  signal-to-noise  ratio  of  the  A-mode 
(attenuated)  echo  could  be  assured  to  be  no  greater  than  unity  (0  db) . 

This  second  threshold  level  would  permit  a 0.66  probability  of  detecting 
echoes  with  an  SNR  of  7 (8.5  db)  and  would  allow  false  alarms  with  a 
probability  of  .002  per  trial,  or  .05  per  search  of  a 25-cell  search 
area.  These  probabilities  are  not  far  from  those  we  initially  assumed 
as  reasonable  for  decisions  involving  a two-of-five  decision  criterion. 

Summarizing  the  above  discussion,  we  have  considered  the  threshold 
levels  required  to  determine  whether  or  not  echoes  are  present  in  the 
vicinity  of  a softly  bounded  MOR.  For  the  example  given,  a two-  out-of- 
five  decision  logic  was  used  to  obtain  decisions  with  a 95%  probability 
of  correctness.  A 4 a threshold  level  was  found  to  be  required  to 
reject  echoes  of  SNR < 3 (4.8  db)  with  the  required  probability  (Pp^  = 0.07) 
and  a threshold  level  of  3.5  a was  found  necessary  for  SNR < 1 (<0  db) . 

These  threshold  levels  provide  a 0.66  probability  of  detection  (as  re- 
quired for  the  95%  probability  of  achieving  two  or  more  detections  ir. 
five  trials)  for  signals  with  SNR's  of  9.7  db  and  8.5  db,  respectively. 

Thu  probability  of  a false  alarm  caused  by  noise  alone  in  a 25-cell 
tracker  acceptance  region  was  found  to  be  0.0075  and  0.05,  respectively. 
Consequently,  the  use  of  the  4 o threshold  permitted  the  critical  region 
of  the  B mode  MOR  to  be  defined  as  that  lying  outside  the  4.8  db  SNR 
contour  of  the  A mode's  MOR,  but  the  use  of  this  threshold  level 
desensitizes  the  radar  by  about  1 db.  If  the  3.5  a threshold  is  used, 
the  critical  region  must  be  limited  to  the  region  outs;de  the  0 db 
contour  of  A mode  propagation,  but  there  would  result  negligible  de- 
sensitization of  the  radar. 

The  above  example  indicates  how,  by  properly  selecting  the  critical 
region  of  the  f-x  plane  for  mode  identification  tests,  satisfactory  tests 
can  be  performed  even  at  soft  boundaries.  The  difficulty  in  the  practical 
implementation  of  such  tests  is  that  of  determining  the  SNR  contours.  Where 

it  is  necessary  to  determine  the  boundary  with  any  degree  of  exactness, 
one  must  be  pessimistic  because  there  are  numerous  contingencies  such  as 
polarization,  which  will  confound  such  predictions.  It  would  be  unwise 


to  attest  such  mode  identification  tests,  for  example,  to  decide 
between  high -ray  and  low-ray  propagation  unless  the  unambiguous  region 
in  which  the  test  is  to  be  made  were  extensive  enough  to  permit  test- 
ing well  beyond  the  estimated  position  of  the  selected  SNR  contour  of  the 
high -ray  MOR.  In  our  earlier  discussions  of  the  extent  of  the  MOR's, 
examples  were  given  in  which  the  extent  of  the  F region  MOR  seemed  to  be 
insufficient  for  this  purpose.  This  was  particularly  the  case  where  an 
underlying  layer  obscured  the  F layer  at  frequencies  only  a little  less 
than  the  MUF  for  any  range. 

3. S Application  of  Tests  to  Fixed- Frequency  Radar 

The  kinds  of  tests  described  in  the  above  sections  are  easily 
applied  to  a frequency -hopping  radar  because  the  propagation  to  any 
given  range  via  a particular  mode  is  usually  strongly  frequency  de- 
pendent. A radar  that  restricts  its  operation  to  a single  frequency 
has  only  the  change  in  range  of  moving  targets  to  use  as  the 
variable  in  analyzing  the  mode  of  propagation.  In  general  for  each 
test  made  by  varying  the  radar's  operating  frequency,  there  is  a 
corresponding  one  that  can  be  made  by  observing  the  echoes  from  a 
target  which  varies  in  range.  Unfortunately,  one  of  the  strongest 
tests  using  frequency  variation- -that  of  testing  at  the  MUF — becomes 
awkward  when  attempted  with  range  variation.  The  passage  of  a target 
into  a mode's  skip  region  corresponds  to  increasing  the  frequency  a bove 
the  MUF  to  the  target's  range,  and  the  strength  of  the  conclusion  that 
can  be  drawn  is  the  same  in  both  cases.  However,  the  test  has  the  un- 
fortunate property  of  providing  the  information  to  determine  the  propa- 
gation mode  at  the  very  instant  that  the  target  may  disappear  from  the 
radar's  field  of  view.  Because  this  is  undesirable,  it  is  preferable 
to  attempt  tests  based  on  the  first  appearance  of  an  echo  at  a boundary. 
Such  tests  are  possible  under  the  following  propagation  and  target  con- 
ditions. 


(1)  Receding  targets  at  skip  distance. 

(2)  Approaching  targets  at  boundary  of  maximum  range  focusing 
(hard  boundary) . 

(3)  Approaching  target  at  undetectable  SNR  contour  (soft  boundary). 


3-20 


In  order  to  make  valid  tests  of  this  kind,  it  is  necessary  that 
the  radar's  detection  process  be  applied  continuously  to  the  critical 
region  (in  this  case,  to  the  critical  radar  delay  intervals).  Indeed, 
if  an  m of  n decision  criterion  is  to  be  applied,  the  testing  must  be 
carried  out  over  a sufficient  interval  of  radar  delay  in  the  critical 
region  that  the  necessary  n observations  can  be  made  of  the  moving 
targets.  A similar  radar  delay  interval  of  detection  processing  can  be 
carried  out  just  within  the  ambiguous  region.  Echoes  first  detected 
in  the  unambiguous  region  can  then  be  tracked  unambiguously  into  the 
region  where  ambiguity  would  otherwise  exist.  If,  on  the  other  hand, 
the  echo  is  detected  within  the  ambiguous  region  without  first  having 
been  detected  in  the  unambiguous  region,  it  can  be  assigned  to  the 
second  mode  (that  at  whose  MOR  boundary  the  testing  is  performed), 
figure  3-7  illustrates  the  procedure. 

If  the  testing  is  to  be  performed  at  a soft  boundary,  considerations 
similar  to  those  discussed  in  the  prior  section  apply.  The  critical  region 
must  be  limited  to  a region  where  the  SNR  of  the  attenuated  mode  is  below 
the  level  used  in  the  design  of  the  decision  process. 

Clearly,  this  kind  of  mode  identification  procedure  will  suffer  if 
track  continuity  cannot  be  maintained  and  so  may  only  be  feasible  with 
trackers  that  operate  with  high  reliability.  This  means  the  tracker  must 
include  adequate  provision  for  carrying  the  track  through  the  fades  that 
commonly  occur  on  HF  signals. 


RADAR  DELAY. 


TIME  — 


Figure  3-7.  Use  of  Range  Variation  by  Target  in  Mode  Identification 
Test.  Approaching  targets  first  appear  at  delay  interval 
T if  B-mode  supports  propagation,  and  at  delay  interval 
t , if  A-mode  supports  propagation. 
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4.  AMPLITUDE  FLUCTUATIONS  AS  AN  INDICATOR 


OF  PROPAGATION  MODE 


In  the  past,  operational  radars  have  used  observed  variations  in 
the  amplitude  of  the  skywave  backscatter  as  a means  for  sorting  out  propa- 
gation modes.  The  information  gained  by  this  process  was  used  for 
propagation  management,  i.e.,  to  insure  that  the  region  of  the  frequency- 
radar  delay  plane  which  was  being  displayed  to  the  operator  mapped  to 
the  geographical  region  of  interest  to  the  users.  This  kind  of  infor- 
mation could  be  extracted  from  the  behavior  of  backscattered  echoes 
from  the  earth's  surface  provided  the  characteristics  of  the  ground 
clutter  were  at  least,  in  part,  determined  by  the  propagation  mode. 

For  the  purposes  of  this  study,  this  kind  of  information  has  been 
assumed  to  be  already  available,  and  we  have  been  concentrating  on  the  i 

question  of  how  to  identify  the  mode  by  which  a particular  target 
echo  has  propagated,  when  the  echo  appears  in  a region  of  the  frequency- 
radar  delay  plane  where  it  is  known  the  echoes  that  have  propagated  via 
two  or  more  modes  can  appear.  It  is  the  purpose  of  this  section  to 
examine  whether  each  propagation  mode  imprints  its  own  character* 
istic  on  the  amplitude  fluctuations  observed  on  an  echo  from  a discrete 
echo  in  the  same  manner  as  has  been  observed  in  the  echoes  from 
azimuthal ly  and  radially  distributed  ground  clutter.  If  such  a de- 
termination can  be  made,  it  would  be  of  great  utility  for  the  identifi- 
cation of  the  mode  operative  in  the  observation  of  echoes  with  a fixed 
frequency  radar,  which  otherwise  only  has  the  ability  to  test  for  mode 
identification  when  the  target  enters  or  leaves  the  geographical 
region  covered  by  a given  mode. 

We  first  discuss  the  amplitude  fluctuation  effects  which,  in  the 
past,  have  been  observed  and  utilized  for  profile  management.  A dis- 
cussion of  -lie  mechanisms  responsible  for  the  fading  of  ground  clutter 
is  next,  followed  by  a discussion  of  how  these  mechanisms  operate  on  a 
discrete  target.  Finally,  a section  discussing  the  observed  fluctua- 
tions of  signals  from  discrete  scatterers  is  presented. 
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4. 1 Use  of  Clutter  Fading  to  Determine  Mode  Coverage 

As  a result  of  an  extensive  series  of  experiments  performed  in 
1965-1966,  it  was  found  that  the  E and  F modes  could  be  distinguished  by 
the  fade  rates  that  appeared  on  the  ground  clutter  echoes.  It  was  found 
that  fixed- frequency  observations  indicated  a fading  period  in  the  ground 
clutter  that  spanned  a few  seconds  (1  to  5 seconds)  when  the  F layer  sup- 
ported propagation  and  some  tens  of  seconds  (15  to  30  seconds)  when  the 
propagation  was  supported  by  the  E layer.  Thus,  an  examination  of  Z- 
axis  records  (echo  strength  displayed  as  recording  density,  as  a 
function  of  radar  delay  and  time)  made  with  slow  frequency  scanning 
(dwells  of  approximately  one  minute  on  each  frequency)  permitted  a 
determination  of  the  regions  of  the  frequency- radar  delay  plane  in 
which  either  E -layer  or  F-layer  propagation  was  dominant. 

This  method,  when  augmented  with  information  extracted  from  an 
examination  of  meteor  echoes  and  of  discrete  ground  targets  of  opportun- 
ity provided  the  means  by  which  radar  profiling  was  accomplished  without 
having  recourse  to  a cooperative  beacon.  Although  the  methods  have 
only  been  applied  successfully  to  the  problem  of  distinguishing  oetweer. 
the  modes  supported  by  the  sporadic  E layer  and  the  F layer,  users 
of  the  method  express  optimism  th3t  it  can  similarly  be  made  to  work 
where  the  mode  observable  regions  of  the  FI  and  F2  layers  are  to  be 
distinguished  and  in  distinguishing  among  other  corntmati  as  of  modes 
as  well. 

• - Fading  Mechanisms 

Fading  i>  a time  variation  in  the  amplitude  of  a received 
signal.  Such  time  variations  can  be  produced  by  two  general  classes 
of  mechanisms.  The  first  is  u reduction  of  the  signal  power  received 
because  of  attenuation  of  the  signal  introduced  on  the  propagation 
path.  The  second  general  class  of  processes  involves  the  interference 
among  wrves  that  have  propagated  over  paths  of  differing  phase 
lengths.  \ CM  signal  will  undergo  this  interference  fading  because  of 
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the  combining  of  signals  that  have  propagated  via  the  differing  propa- 
modes.  However,  here  we  presume  such  modes  to  have  been  time- 
, so  this  kind  of  multi-mode  interference  does  not  affect  the 
amplitude  of  concern.  However,  there  are  several  forms  of  multi- 
path  propagation  which  an  QTH  radar  of  moderate  bandwidth  ( <0.S  MHz)  does 
not  resolve.  Among  these  are  the  magneto-ionic  mult’moding  and  the  so-called 
fine-grain  multipath  which  results  from  ionospheric  irregularities. 

Also  included  here,  under  the  category  of  interference  effects, 
will  be  those  caused  by  the  interference  of  multiple,  but  unresolved, 
scatterers. 

4.2.1  Fading  Caused  by  Signal  Attenuation  (Energy  Loss) 

There  are  three  mechanisms  by  which  an  ionospherically  propa- 
gated signal  can  become  significantly  attenuated;  these  are 
absorption  of  the  signal,  defocussing  of  the  signal,  and  scattering 
of  the  signal  by  interposed  irregularities. 

Absorption  occurs  predominantly  in  the  D layer,  through  which 
HF  skywave  signals  all  must  pass,  regardless  of  their  mode  of  propa- 
gation. This,  along  with  the  fact  that  the  fluctuations  in  absorption 
occur  rarely  enough  that  their  availability  cannot  be  counted  on  for 
mode  identifying  tests,  would  seem  to  rule  out  their  use.  This  is 
particularly  the  case  for  the  sudden  ionospheric  disturbances  (SID's) 
which  are  produced  throughout  t’.„.;  sunlit  side  of  the  earth  by  solar- 
flare -associated  x-rays. 

Ionization  in  the  absorptive  layers  can  also  be  produced  by 
energetic  particles,  either  in  association  with  auroral  phenomena 
or  within  the  po  ar  cap  where  low  energy  solar  protons  have  access  to 
the  absorptive  layer.  Because  of  the  spatial  localization  of  the 
ionization  so  produced,  the  possibility  exists  that  some  kind  of  mode 
identification  could  be  based  upon  complete  knowledge  of  the  spatial 
distribution  of  the  absorption  and  of  the  ray  paths  of  the  propagation 
nodes.  However,  we  are  a long  way  from  possessing  sufficiently  detailed 
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information  regarding  the  distribution  of  auroral  ionization  to  make  this 
approach  feasible. 

Defocus ing  of  the  signal  can  occur  either  because  of  a 
vertical  distribution  of  an  otherwise  uniform  electron  density  (as 
discussed  in  Section  2)  or  as  the  result  of  some  horizontally 
structured  perturbation  of  electron  density  such  as  can  be  produced  by 
internal  gravity  waves  in  the  atmosphere.  The  former  mechanism  would 
produce  soft  boundaries  to  mode-observable  regions,  as  discussed 
earlier. 

Internal  gravity  waves  in  the  atmosphere  can  also  produce  electron 
density  inhomogeneities  with  horizontal  and  vertical  structure,  which 
focus  and  defocus  the  radar's  illumination.  These  effects  will  have 
characteristic  periods  of  about  fifteen  minutes  or  longer.  Irregulari- 
ties . ssociated  with  air  motions  of  shorter  periods  have  also  been 
observed.  Waves  with  periods  less  than  five  minutes  are  refern d to 
as  buoyancy  waves  cr  acoustic  waves,  the  distinction  being  based  on  the 
relative  magnitudes  of  various  terras  in  the  equation  of  motion.  For 
our  present  purpose,  we  only  note  that  these  perturbations  of  electron 
density  occur  in  the  1:  and  F layers  and  that  the  attendant  focussing 
and  defocuss’ng  is  one  of  several  causes  of  fluctuation  to  be  dis- 
cussed subsequently. 

Signal  strength  can  also  be  decreased  by  scattering.  One 
phenomenon  of  this  category  is  the  obscuration  by  the  sporadic  H layer. 
Another  is  the  similar  obscuratim  caused  by  auroral  ionization.  How- 
ever, when  such  an  effect  exists,  a more  noticeable  signal  fluctuation 
occurs  through  interference  phenomena,  to  be  discussed  in  a later 
section.  Again,  we  note  only  that  scattering  is  one  of  the  phenomena 
that  contribute  to  the  fluctuation  in  signal  strength. 

4.2.2  Amplitude  Fluctuation  by  Interference  F.ffects 

Interference  effects  are  responsible  for  much  of  the  signal 
strength  fluctuation  observed  in  the  skywavc  signals.  An  0 11 1 radar  of 
moderate  bandwidth  (20  kHz)  is  immune  to  some  of  the  effects  encountered 


in  other  services  utilizing  narrow  bandwidths  because  the  signals  that 
have  propagated  via  the  gross  modes  (e.g.,  E layer,  FI  layer,  F layer 
high- ray,  F layer  low -ray)  are  time  resolved  and  so  do  not  interfere 
with  one  another.  However,  finer  mode  structure  such  as  the  magneto- 
ionic raultimoding  is  often  not  resolved  and  so  will  manifest  itself 
as  an  interference -caused  fluctuation  of  signal  strength.  Irregulari- 
ties are  also  a source  of  fine  grain  multipath,  but  their  effect  is 
more  frequently  analyzed  in  terras  of  a random  diffracting  screen 
(e.g.,  Booker,  Ratcliffe,  and  Shinn,  1949).  Extensive  literature 
exists  on  this  topic  in  which  the  fading  is  referred  to  as  "scintilla 
tion." 

The  magneto- ionic  multimoding  occurs  most  strongly  in  signals 
that  have  propagated  via  the  F- layer  where  the  phase  path  through  the 
magneto- ionic  medium  is  very  long.  The  irregularities  have  been 
observed  at  both  E layer  and  F layer  heights.  The  F layer  irregulari- 
ties would,  of  course,  affect  only  F layer  modes,  while  the  E layer 
irregularities  could  introduce  fluctuations  on  all  modes. 

Time  variations  are  required  in  each  case.  The  magneto-ionic 
multimodes  will  change  their  relative  phase  as  the  layer  shape,  and 
consequently  the  ray  oath,  geometry  varies.  Wind-induced  drifts  are 
usually  invoked  to  explain  the  time  variation  of  the  irregularities 
in  the  E layer  anc  magnet cspheric  electric  fields  to  explain  F -layer 
irregularity  drifts. 


These  fluctuations  may  be  regarded  as  having  their  origin 
entirely  in  the  ionosphere.  That  is  to  say,  they  wo<Jld  occur  regard- 
less of  the  nature  of  the  scattering  target,  be  it  a discrete  or  a 
distributed  scatterer . We  next  take  up  some  mechanisms  for  signal 
strength  fluctuations  which  depend  upon  the  nature  of  the  target. 


The  first  invokes  time  vari''ions  in  the  properties  of  the 
target  to  account  for  signal  itrength  fluctuations.  A target  of 
sufficient  size  (several  wavelengths)  will  have  a scattering 
cr os/  section  that  varies  with  the  aspect  from  which  it  is  viewed  by  a 
radar.  Multiple,  unresolved  targets  will  scatter  hack  echoes  tha  «y 


interfere  constructively  or  destructively,  and  so  changes  in  their 
relative  position  will  produce  signal  fluctuations.  A similar  effect 
occurs  between  a moving  target  and  its  image  in  the  ground.  Also, 
sea  backscatter  has  a fluctuation  as  the  approaching  ana  receding 
Bragg-scattered  components  interfere,  and  as  higher  order  scatterers 
interfere  with  each  other  and  the  first  order  components.  The  above 
mechanisms  are  all  similar  in  requiting  a motion  of  the  scatterer 
to  operate. 


Consider  next  the  backscatter  from  a distributed,  but  fixed, 
target  such  as  from  ground  clutter.  A pulsed  radar  with  a time  resolution 
leading  to  an  interval  of  confusion  At,  and  having  azimuthal  resolution  A0, 
responds  to  surface  features  by  an  amount  determined  by  performing  a 
vector  integration  over  all  scatterers  in  the  illuminated  field. 


Thus,  the  area  including  the  scatterers  is 


A (r,  <$>)  = r A6  c At  seed 


N,  the  number  of  these  scatterers  is  proportional  to  A.  The  ,summat ion 
over  the  set  of  scatterers  in  an  unresolved  area  it  mathematically 
equivalent  to  an  N’-step,  two-dimensional  random  walk  with  the  consecpience 
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that  the  resulting  amplitude  has  :ero  mean  and  a variance  proportional 
to  N.  This  variance  is  observed  as  the  returned  power.  Thus,  the  use 
of  a ground  clutter  coefficient  5 with  the  illuminated  area  to 
determine  the  clutter's  radar  cr-'ss  section  is  justified,  and  we  can 
write  the  returned  power  or  variance  of  the  received  amplitude  as 

P “ o c A r Ad  Ax  sec  « 

A further  character icat ion  of  the  clutter  field  is  required 
for  our  purposes.  This  is  an  indication  of  the  sensitivity  of  the 
clutter  returned  from  any  given  range  to  changes  in  the  angle  of 
arrival  of  the  radar  illumination.  Note  that  we  arc  not  interested 
here  in  the  sec  factor  in  the  expression  for  the  returned  clutter 
power.  Rather,  wc  are  Interested  in  how  quickly  the  random  amplitude 
can  change  from  a maximum  value  to  a low  value  as  the  angle  a is 
changed.  This  characteristic  is  measured  by  the  autocovariance  of 
the  angular  spectrum  of  the  clutter  returns  from  a resolved  area,  lor 
our  purpose  we  may  obtain  an  estimate  of  this  angular  sensitivity  by 
identifying  the  quantity  with  the  resolution  of  an  nr  ay  with  an 
iperturc  equal  to  the  depth  of  the  unresolved  cell,  1)  = e.\i  secu.  The 
effective  wavelength  is  half  the  electromagnetic  wavelength  because  any 
1 1 -placement  along  the  aperture  affects  the  phase  of  the  incident  and 
b.n  kscattered  wa».e<. 
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Iho  difference  in  phase  between  the  signal  of  the  front  half  of 
the  cell  to  that  of  the  back  half  is 


’iC  At  seca 
A<t>  = 2n  *i  A 

= 2tt  f At  seca 


Taking  the  derivative,  wc  have 


dA4> 

da 


= 2n  f Ai  seca  tana 


And  solving  for  the  increment  of  a that  will  produce  a phase  change  of  it  radians 
we  obtain 


An  • 

v da 


1 


2 f At  seca  tana 


cosa  cota  _ cos  a 


2 f At 


2 fATsina 


Values  of  Aa  are  plotted  for  various  frequencies  and  values  of 
a in  Table  4-1.  A radar  resolution,  At,  of  50  ps  is  assumed. 


Table  4-1.  Width  of  Angular  Covariance  (degrees) 


<_(Mllz) 

5 

7 

10 

15 

20 

25 

30 

5 

1.7 

0.93 

0.65 

0.43 

0.52 

0.26 

0.21 

10 

0.64 

0.45 

0.32 

0.21 

0.16 

0.12 

0.10 

15 

0.41 

0.29 

0.20 

0.15 

0.10 

0.08 

0.06 

20 

0.30 

0.21 

0.15 

0.10 

0.07 

0.06 

0.05 

25 

0.22 

0. 16 

0.11 

0.07 

0.06 

0.04 

0 . 04 

50 

0. 17 

0.  12 

0.09 

0.06 

0.04 

0.03 

0 . 03 

55 

0. 13 

0. 10 

0.07 

0.04 

0 03 

0.02 

0.02 

10 

0. 10 

0.07 

0.05 

0.03 

0.03 

0.02 

0.02 

45 

0.08 

0.06 

0 . 04 

0.03 

0.02 

0.02 

0.01 

50 

0.0b 

0.04 

0.03 

0.02 

0.02 

0.01 

0.01 
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It  is  clear  from  Table  4-1  that  rather  small  changer,  in  the  angle  at 
which  a resolvable  cell  is  observed  with  the  radar  can  produce  large 
changes  in  the  level  of  the  radar  echo.  Also,  note  that  the  sensitivity 
to  the  angle  of  arrival  decreases  by  more  than  an  order  of  magnitude  over  the 
interval  5°  to  50°.  The  implication  of  this  is  that,  for  a given  rate 
of  ionospheric  tilting,  such  an  effect  would  lead  to  more  rapid 
fluctuations  of  the  ground  clutter  observed  at  near  ranges  than  those 
observed  at  more  distant  ranges,  Mitra  (1949)  presents  some  curves  of  the 
virtual  height  of  the  F,  FI,  and  F2  layers,  as  functions  of  the  time 
of  day,  from  which  can  be  scaled  an  estimate  of  the  ionospheric  curva- 
ture (change  of  layer  tilt)  in  the  east-west  direction. 

2 

Let  h=a  + b,  x + c.x  - 
o 1 1 

= b + 2c,  x 
dx  1 1 

d2h  _ - 
~ 2cl 

where  x is  east-west  distance  and  where  aQ,  b^  and  Cj  are  appropriate  expan- 
sion coefficients  to  fit  the  curve  for  Washington.  D.C.  summer  presented  by 
Mitra  (1952). 

Near  sunrise  the  layer  height  descends  at  a rate  of  SO  km/hour. 

About  two  hours  later,  this  rate  has  decreased  to  no  more  than  10  km/ 
hour.  Thus,  we  have 


d“h  „ 40  km/hour 

dx?  " 2 hours 


20  km 
hour-' 


At  the  latitude  of  Washington,  an  hour  corresponds  to  a distance 

x = — • 2n  r • cos  L = 1348  km. 

24  e 

where  L is  the  latitude,  and  rg  is  the  radium  of  the  earth.  Wc  may  change 
the  units  of  the  above  estimates  of  curvature  to 


d~h  .radians. 
dx2  *■  sec. 


20  km  # 1 hour  _ hour 

hour^  1348  km  3600  sec. 
, , .-6  radians 


Or  in  units  of  degrees,  about  2.30  x 10  degrees/second.  Were  this 
the  only  component  of  layer  curvature,  the  characteristic  time  for  a 
fade  produced  by  this  mechanism  would  be  greater  than  40  seconds. 


0.01  degrees 
2 . 3x10" 4 degTsec. 


> 40  seconds 


The  gravity  waves  discussed  earlier  impose  additional  components 
of  ionospheric  structure  that  can  lead  to  variations  in  the  angle  of 
arrival  of  the  radar  illumination  on  target.  Clark,  et  al.  (1971)  have 
performed  modei  calculations  of  the  perturbation  of  the  F layer  electron 
density  produced  by  gravity  waves  of  various  periods  and  wavelengths. 

As  one  example,  the  effects  of  an  eastward  propagating  wave  with  a 
wavelength  of  about  300  km  (k  - 0.02  km  *)  and  period  of  20  minutes 
was  computed.  The  amplitude  of  the  wave  was  2 x 10  dyn*  sec  , which 
led  to  a wavelike  deformation  of  the  contours  of  equal  electron 
density  with  peak-to-peak  amplitudes  of  order  50  km.  Idealizing  the 
effect  as  that  of  producing  a reflecting  surface  that  deviates  from 
a height  of  300  kra  according  to  expression 


h = A cos  k x 

permits  us  to  determine  the  shape  of  the  layer,  ^ . Thus, 


-j—  = - A k sin  k x 

dx 

Then  taking  a time  derivative  produces  the  rate  of  change  of  slope. 


_d  rdh'| 
dt  (dxj 


. ,2  , dx 

- A k cos  k x T- 
at 


dx 

The  gravity  wave's  velocity  is  equal  to  its  length,  A, 
divided  by  its  period  T.  In  this  example,  A = 300  km,  and  T = 20  min. 
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therefore 


dx 

dt 


A 300  km 
T ”20  mins. 


= 15  km/minute 


The  maximum  magnitude  of  the  rate  of  change  of  slope  occurs 
when  we  have 


I—  f— 1 
'dt  IdxJ 


dhl | = A k 


2 dx_ 
dt 


max 


= 25  km  x (0.02) 2 km"2  15  ^ , 

minute 

= 0.15  radians/rainute 

_ 3 

= 2.5  x 10  radians /second 

= 0.143  degrees/second 

Referring  to  Table  4-1,  we  see  that  except  at  the  lowest  frequencies 

t=S  MHz)  incident  at  the  lowest  elevation  angles  (-5°),  an  angular 

drift  of  this  magnitude  is  sufficient  to  produce  fading  with  periods 

less  than  a second.  Clark,  et  al.  (1971)  computed  the  electron 

density  perturbation  produced  by  a model  gravity  wave  with  length 

about  1600  km  and  period  60  minutes.  The  amplitude  of  the  gravity  wave 

9 U 2 

in  this  case  was  2 x 10  dyr. ‘ sec  . This  wave  produced  a corrugation 
of  the  constant  electron  density  contours  of  about  20  km  peak  to  peak 
amplitude  at  the  F layer.  These  values  lead  to  a maximum  rate  of  angular 
deviation 


\± 

' dt 


dh) 

dx 


max 


i A k2  — 
A K dt 


= 10  k.n  (.004)2  km  2 * km/minute 

uU 

= 0.00426  radians/minute 
= 0.004  degrees/second 

Again,  referring  to  Table  4-1  we  find  that  these  waves  could  produce 
fading  with  a time  scale  measure  in  seconds  or  tens  of  seconds  over 
much  of  the  frequency  and  elevation  angle  interval  of  interest. 


4.3  Fading  Mechanisms  as  They  Affect  Echoes  from  Discrete  Targets 

The  fading  of  echoes  from  a discrete  target  are  of  interest  to 
us  for  two  reasons.  First,  a tracker  operating  on  an  OTH  radar  must  be 
able  to  bridge  the  periods  of  low  detection  probability  that  occur 
during  fades.  Secondly,  we  wish  to  determine  whether  the  fading  rate 
of  the  echo  is  characteristic  of  the  propagation  mode  in  the  same 
manner  as  has  been  observed  in  the  fading  of  ground  clutter. 

Not  all  the  fading  mechanisms  discussed  in  the  preceding 
section  would  contribute  to  the  fading  of  a discrete  target.  Particu- 
larly to  be  excluded  in  this  respect  would  be  those  that  depend  upon 
the  target  being  distributed,  i.e.,  those  involving  a change  in  the 
relative  position  of  unresolved  scatterers.  Of  those  that  do  affect 
echoes  from  discrete  targets,  not  all  are  capable  of  providing  mode 
identifying  information.  Of  this  kind  would  be  the  fades  caused  by 
sudden  ionospheric  disturbances,  which  similarly  affect  all  modes. 

Some  of  the  mechanisms  are  of  a kind  chat  would  differently 
affect  E-mode  and  F-mode  signals.  Magneto-ionic  multimoding  is  one  of 
these.  For  the  E mode  the  change  with  range  of  the  difference  in  phase 
between  the  ordinary  and  extraordinary  modes  is  rather  small,  as  mani- 
fested by  the  small  Faraday  rotation  observed  in  E-mode  signals 
(Earhart,  1973).  That  for  the  F mode,  on  the  other  hand,  is  large. 
Consequently,  were  this  the  sole  source  of  fading  in  discrete  target 
echoes,  mode  identification  for  major  modes  would  be  simple. 

The  effect  of  gravity  waves  is  far  more  pronounced  in  the  F 
layer  than  in  the  r.  layer;  so  again,  we  can  anticipate  a distinct 
difference  according  to  propagation  mode  from  this  mechanism  as  well. 

The  remaining  mechanisms  are  those  associated  with  ionospheric 
irregularities.  In  most  analyses,  the  effects  of  irregularities  are 
regarded  as  random,  as  opposed  to  being  analyzable  into  regular  wavelike 
or  periodic  structures  used  to  describe  the  magneto-ionic  multimode  inter 
fercnce  (Faraday  rotation)  and  the  focusing  effects  of  traveling  ionosphe 
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disturbances.  Also,  the  structure  of  the  effects  is  rather  different.  While 
irregularities  structure  the  radar's  illumination  fields  into  interference- 
caused  "speckles,"  measuring  a kilometer  or  less,  the  other  mechanisms  cause 
focusing  or  polarization  bands  that  measure  tens  or  hundreds  of  kilometers 
in  width.  If  the  effects  are  superposed,  the  irregularity  fading  is  apt  to 
mask  the  other  effects. 


The  theory  of  ionospheric  irregularities  that  cause  signal  scintil- 
lation is  not  as  well  developed  as  those  for  Faraday  rotation  and  traveling 
ionospheric  disturbances.  Particularly  lacking  are  theories  which  can  predict 
differences  in  the  scale  size  or  probability  of  occurrence  that  may  exist 
between  E- layer  and  F-layer  irregularities.  For  this  reason,  the  following  sec 
tion  examines  empirically  the  difference  in  fading  rates  that  occur  in  E-and 
F-layer  modes. 


4 . 4 Empirical  Evidence  Regarding  Usefulness  of  Target  Fading 

in  Mode  Identification 

Although  abundant  literature  exists  on  the  topic 
of  fading  HF  signals,  the  present  study  must  be  restricted  to  those 
effects  which  persist  when  the  time  resolution  of  the  probirg  signal  is 
sufficient  to  resolve  the  multiple  propagation  modes.  That  is  to  say, 
we  must  exclude  from  consideration  those  effects  which  are  themselves 
manifestations  of  multimode  propagation  that  would  be  resolved  by  the 
radar . 

Experimental  evidence  has  been  collected  by  sounding  on 

closely  spaced  frequencies,  indicating  that  signal  fading  is  associated 

with  effects  that  occur  near  the  reflection  height.  (Essex  and  Hibberd, 

1968,  and  Essex,  1968).  Thus,  if  the  physical  processes  responsible  for 

the  observed  echo  fading  were  significantly  different  for  the  different 

« regions  of  the  ionosphere,  it  would  be  reasonable  to  expect  that  each 

* . 

j region  imparts  a characteristic  signature  or  modulation  on  a signal  which 

could  be  used  for  made  identification. 

However,  the  use  of  the  fading  characteristic  of  the  target 
for  mode  identification  is  complicated  by  the  fact  that  the  fading 
is  usually  the  result  of  time  variations  of  a spatially  structured 
illumination  pattern.  Thus,  motion  of  the  target  would  modify  the 
1 effects  that  have  been  observed  over  oropagation  paths  with  fixed 

t 

end  points.  If  the  target  were  motionless  or  moved  sufficiently  slowly, 
the  characteristics  observed  by  fixed  end-point  sounding  would  be 
; immediately  applicable.  On  the  other  hand,  if  the  targets  moved  suf- 

ficiently rapidly,  the  expected  fading  characteristics  could  be  pre- 
dicted by  considering  the  rate  at  which  the  target  passes  through  an 
essentially  stationary  or  "frozen"  illumination  pattern.  Considerably 
greater  difficulty  is  to  be  expected  when  the  fade  rates  associated 
with  changes  in  the  propagation  path  are  comparable  to  those  caused 
by  the  motion  of  the  target  through  the  spatially  structured  illumina- 
tion. As  we  shall  see,  this  latter  case  is  all  too  often  encountered. 
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In  this  section  the  means  and  requirements  for  performing 
mode  identification  are  discussed.  The  basis  of  this  discussion  is 
this  Laboratory's  experience  and  the  information  available  from 
published  literature.  Many  mechanisms  have  been  proposed  to  describe 
and  characterize  the  observed  fading  effects.  Some  of  these  mechanisms 
have  been  discussed  previously.  The  purpose  of  this  section  is  not  to 
validate  the  theories  ani  mechanisms  invoked  to  explain  the  observa- 
tions of  signal  fading;  instead  the  purpose  is  to  discuss  and  indicate 
possible  ways  in  which  the  fading  of  the  signal  might  be  applied  to 
solve  the  problem  of  mode  identification. 

4.4.1  Background  and  Measurement  Techniques  Used  to  Characterize 

Amplitude  Fading  of  Signals 

The  temporal  variation  of  signals  can  be  determined  by  measur- 
ing the  received  amplitude  of  each  pulse  transmitted  over  the  path.  The 
temporal  variation  or  fading  of  ionospherically  propagated  signals  has 
been  investigated  extensively,  dating  back  to  Ratcliffe  and  Pawsey 
(1933) . Considered  here  are  only  those  observations  made  with  pulse 
length  short  enough  (less  than  200  usee.)  that  the  transmissions  via 
the  different  modes  are  resolved.  The  ordinary  and  extraordinary 
components  can  be  separated  on  the  Hsis  of  time-of-arrival  or  delay 
uime  with  sounders  having  time  resolutions  of  one  microsecond  order. 
Additionally,  a number  of  investigators  utilized  circularly  polarized 
antennas  to  further  ensure  that  the  cause  of  the  fading  was  not  Faraday 
rotation. 

Although  fading  produced  by  Faraday  rotation  could  be  used 
as  an  indicator  of  propagation  mode,  this  section  stresses  fading 
caused  by  small  scale  ionospheric  irregularities  which  are  changing  with 
time  and  might  be  different  for  the  different  modes.  These  effects 
occur  frequently  enough  ana  obscure  the  more  predictable  Faraday 
rotation  so  that  they  must  be  contended  with  in  the  development  of  a 
mode  identification  procedure.  As  indicated,  some  of  the  other  causes 
of  the  signal  fading  can  be  eliminated  by  proper  choice  of  operating 
parameters  so  that  the  cause  of  fading  can  be  limited  largely  to  tnoso 
cases  of  interest. 


W 
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Fading  effects  have  been  characterized  and  classified  in 
various  ways.  Of  particular  interest  to  us  are: 

o Fading  frequency  spectrum  or  fade  rate. 

o Amplitude  probability  density  function 

Each  of  these  methods  has  advantages  and  offers  a different 
dimension  for  studying  the  amplitude  fading  of  a signal.  The  spectrum 
analysis  separates  the  fading  rate  into  the  various  fading  frequency 
components.  Sometimes  a crude  estimate  of  the  fading  spectrum  is  made 
by  counting  the  number  of  observed  amplitude  maxima  (or  minima)  per 
unit  of  time. 

The  amplitude  probability  density  function  is  typically 
studied  to  identify  which  type  of  distribution  function,  such  as  Rayleigh, 
Rice,  or  Gaussian,  characterizes  a set  of  measured  amplitudes.  The 
form  of  the  distribution  function,  besides  being  of  interest  in  coamuni ca- 
tions system  design,  implies  certain  physical  characteristics  about  the 
properties  of  the  transmission  medium  such  as  the  strength  and  multi- 
plicity of  the  scattered  components. 


4.4.2 


Previous  Work 


The  reported  data  on  signal  fading  exists  for  two  different 
sounding  techniques,  one  using  essentially  vertical  incidence  propa- 
gation, and  the  other  oblique  incidence  propagation.  Evidence  (McNicol , 
R. , 1949,  and  Kokourov,  et  al.,  1971)  indicates  that  the  effects 
measured  on  oblique  paths  are  different  from  those  measured  simultan- 
eously on  a vertical  path.  Of  interest  for  the  current  application 
to  an  OTH  radar  is,  of  course,  an  oblique  path.  Unfortunately,  the 
larger  quantity  of  data  available  is  for  the  vertical  incidence  propa- 
gation case. 


The  large  volume  of  vertical  incidence  data  was  collected  to 
measure  the  motion  of  ionospheric  winds  and  the  structure  of  the  wind- 
driven  moving  irregularities.  These  investigations  usually  follow  the 
simple  method  of  spaced  receiver  by  Mitra  (1949)  or  the  more 
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sophisticated  techniques  applied  by  Booker,  Ratcliffe,  and  Shinn  (1950) 
or  Briggs,  Phillips,  and  Shinn  (1950).  Although  not  directly  applicable, 
it  is  possible  to  gain  insight  from  the  extensive  amount  of  data  collected 
using  vertical  soundings. 

The  case  of  vertical  incidence  propagation  is  discussed  first 
in  the  following  section.  The  next  section  discusses  the  data  avail- 
able on  fading  for  the  oblique  propagation  case. 

4.4. 2.1  Fading  Characteristics  of  Signals  Observed  at  Vertical 

Inc idence 

Amplitude  distributions  were  measured  as  a function  of  mode 
on  frequencies  of  2 and  4 MHz  by  Flood  (1954).  The  amplitude  fading 
was  categorized  for  each  observation  as  being  either  a Rayleigh  or  a 
Gaussian  distribution.  The  data  were  then  tabulated  and  the  percentage 
of  total  observations  having  each  type  of  distribution  was  reported 
as  a function  of  propagation  mode  as  shown  below: 


r— ""  — ■ ■—  - — 
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l PERCENTAGE  OF 

OBSERVATIONS 

WITH 

MODE 

RAYLEIGH  DISTRIBUTION 

GAUSSIAN  DISTRIBUTION 

H 

59% 

41% 

F 

o\o 

CO 

rr 

52% 

Although  there  is  some  difference  in  the  predominant  type  of  distribution 
for  the  two  modes,  it  can  be  seen  from  this  data  that  the  differences 
are  too  small  to  provide  a means  of  mode  identification. 

Data  on  the  fade  rates  is  often  not  reported  directly,  although  Mich 
data  have  been  used  to  deduce  the  results  indicated.  The  fade  rate  can 
be  constructed  from  quantities  reported  in  the  literature  by  c 'bining 
the  reported  scale  size*and  drift  velocities  of  the  irregularities. 

A large  volume  of  data  has  been  reported  on  rhe  measured 
drift  velocity  of  the  ionospheric  irregularities.  The  measured  velocity 
of  the  winds  varies  with  a number  of  geophysical  factors  including  season 


•:  *an  estimate  o?  mean  • <~.c 
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and  time  of  day.  The  reported  wind  speed  measurements  generally  fall 
in  the  range  of  40  to  150  in/ sec.  Davies  (1969)  quotes  the  average  com- 
ponent of  the  wind  speed  as  being  80  m/sec  in  the  E region  and  about  100 
m/sec  in  the  I--  region.  Rao  and  Rao  (1957)  measured  a wind  speed  of  R4m./sec 
in  the  E region  in  agreement  with  the  reported  results  of  Davies. 

Mitra  (1949)  reports  the  most  frequent  wind  velocity  observed  was 
50m/sec.  Mitra  does  not  segregate  the  data  by  ionospheric  propagation 
mode.  However,  he  specifically  conducted  a series  of  tests  designed 
to  compare  the  wind  speed  in  the  E region  with  that  in  the  F region.  This 
was  done  by  collecting  alternate  samples  from  each  region  using  an 
alternating  delay  time  gate.  On  the  basis  of  a large  number  of  these 
samples,  Mitra  concluded  that  the  magnitude  and  also  the  direction  of 
"wind  was  always  nearly  the  same  for  the  two  echoes"  from  the  different 
regions. 

A more  recent  measurement  of  the  ionospheric  drifts  was 
reported  by  Kokourov  et  al.  (1971);  some  of  their  histograms  are  shown  in 
Figure  4-1.  The  data  are  shown  for  upper  ionosphere  (F  region)  and 
lower  ionosphere  as  measured  at  two  locations  in  the  Soviet  Union 
as  a function  of  season.  The  most  frequently  observed  drift  velocities, 
between  40  to  120  m/sec.  are  essentially  in  agreement  with  previous 
reported  data.  The  important  point  to  be  made  is  that  again  there  is 
no  large  difference  in  the  observed  drift  speed  between  the  E and  F 
regions  of  the  ionopshere.  This  implies  that  the  fade  rate  will  be 
different  in  the  two  regions  only  if  the  scale  sizes  of  the  irregulari- 
ties are  significantly  different. 

Histograms  of  the  frequency  of  occurence  of  different  values 
of  the  scale  size  of  the  irregularities  are  shown  in  Figure  4-2  which  is 
taken  from  Essex  and  Hibherd  (1968).  This  data  indicates  that  the  mean 
size  of  the  E region  and  F region  irregularities  differ  by  a factor  of 
2 to  l.  The  measured  scale  size  of  the  sporadic  E irregularities  is 
between  the  values  measured  for  the  F region  and  E region.  Also,  some 
additional  da:a  were  collected  on  the  Buckland  Park  array  in  Australia 
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Figure  4-2.  Histograms  of  the  Frequency  of  Occurrence  of  Different  Values 

of  the  Horizontal  Scale  Size,  i,  of  the  Ionospheric  Irregularities 


Reproduced  from  Essex  and  Hibberd  (1968) 


•1-20 


by  Felgate  and  Golley  (1971).  The  array  consists  of  80  independent 
aerials.  Tne  output  display  consists  of  light  pattern  formed  by  arraying 
a series  of  lights  in  the  same  manner  as  are  the  antennas  and  driving 
the  output  of  each  light  in  proportion  to  the  signal  strength  observed 
in  the  corresponding  element  of  the  array.  The  observed  light  patterns 
were  grouped  into  seven  classes  that  included  large-scale,  fringes,  crossing 
fringes,  random,  and  no  fading.  The  percentage  of  time  that  each  pat- 
tern was  observed  was  reported.  A comparison  of  these  results  for  the 
large-scale  patterns  and  the  case  of  no  fading  is  shown  for  each  mode 
in  the  table  below. 


rge  Scale 
Fading 


PERCENTAGE  OF  OBSERVATIONS 
FOR  THE 

INDICATED  PATTERN  TYPES 


E 

Es 

12 

19 

2 

4 

The  table  indicates  that  the  fading  of  the  signal  reflected  from  the 
F layer  was  categorized  as  large-scale  or  non-existent  70%  of  the  time 
while  that  from  the  E and  Es  layers  was  so  categorized  only  10  to  20% 
of  the  time.  By  large-scale  the  authors  meant  fading  structure  that  had 
a scale  larger  than  their  1 km  antenna  array  diameter.  Thus,  the  spatial 
structure  of  the  illumination  was  observed  to  have  a scale  less  than 
one  kilometer  in  size  30%  of  the  time  when  the  F layer  produced  the 
reflection  and  in  excess  of  80%  of  the  time  when  the  E layer  produced 
the  reflection.  However,  the  size  scales  for  the  two  modes  do  overlap, 
and  inasmuch  as  the  mode  identification  requires  a differcice  discern- 
ible in  a single  sample  (i.e.,  the  fading  on  a single  track)  rather  than 
being  discernible  in  a statistical  analysis,  these  results  are  not 
conclusively  favorable. 


4.4. 2. 2 Fading  Characteristics  of  Signals  Observed  at  Oblique 

Incidence 

As  previously  indicated,  there  is  less  data  available  on 
fading  statistics  for  the  oblique  propagation  case.  There  are,  however, 
some  limited  results  available  from  the  POLAR  CAP  III  program  and  the 
POLAR  FOX  I program  that,  although  not  extensive,  are  valuable  because 
they  are  directly  applicable  to  the  problem  in  hand. 

In  the  POLAR  CAP  III  program,  a series  of  relative  power 
measurements  were  made  of  echoes  from  a remote  beacon  located  about 
1600  km  from  the  radar.  The  final  report  (ITT-EPL  Staff,  1974). 
discusses  the  results  of  two  case  studies:  one  when  the  propagation 

was  by  a sporadic  E layer  and  the  other  when  the  propagation  was  by  an 
F layer.  The  results  are  shown  for  a 32-pulse  dwell  (one  second)  for 
both  cases  in  Figure  4-3  and  4-4.  The  outputs  shown  are  fully  processed 
data  which  included  correlation  and  Fourier  transformation  for  doppler 
detection.  Each  output  is  an  independent  sample  taken  at  the  doppler 
offset  frequency  of  the  target  as  observed  on  an  amplitude  versus  dop- 
pler frequency  display.  The  sampled  output  was  collected  at  a point 
in  the  processing  chain  at  which  inputs  for  a mode  identification  al- 
gorithm would  be  supplied. 

A total  of  16  seconds  of  data  are  shown  in  Figures  3 and  4. 

The  important  feature  to  note  is  the  large  variation  (18  db)  of  the 
relative  amplitudes  for  the  F mode  case  in  comparison  to  that  (6  db) 
for  the  E mode  case.  It  should  be  emphasized  that  these  data  were 
collected  in  the  polar  environment  and  represents  a very  limited 
quantity  of  data.  There  is  no  means  to  interpret  statistically  how 
applicable  this  data  would  be  to  the  non-polar  region.  This  data 
merely  illustrates  the  type  of  amplitude  variation  differences  nec- 
essary to  perform  mode  identification. 





In  the  data  discussed  up  to  this  point,  the  terminals  used 
for  reception  and  transmission  have  been  located  at  fixed  position. 

Also,  in  the  preceding  examples  the  beacon  was  at  fixed  land  position. 

A typical  airspeed,  say  for  a Boeing  707  aircraft,  is  around  400  mph 
(050  kmh)  or  177  msec.  The  speed  of  the  aircraft  is  only  slightly 
larger  than  the  previously  indicated  measurements  of  the  wind  speed 
within  the  ionopsherc.  Thus,  the  fading  characteristics  of  echoes  from 
an  aircraft  target  could  be  influenced  significantly  by  the  motion  of 
the  aircraft. 

A comparison  of  the  fading  statistics  of  aircraft  and  repeater 
targets  was  made  by  Riverside  Research  Institute  (RRI)  and  was  reported 
at  the  1972  OHD  Technical  Review  Meeting  (Mastrandrea  and  Healy,1972). 
These  data  were  collected  with  a Whitehouse  radar  during  the  POLAR  FOX  1 
experiments.  The  reported  data  show  amplitude  as  a function  of  time  and 
cumulative  distributions  of  the  amplitudes.  The  beacon  was  located 
at  a nominal  range  of  1930  km  from  the  radar.  Two  aircraft  targets 
were  observed  along  the  heading  of  the  beacon  for  the  cases  shown.  The 
observed  shape  of  the  distribution  was  only  slightly  different  (largest 
reported  difference  was  20%)  for  the  combined  results  of  the  two  air- 
craft and  the  beacon  targets.  This  would  appear  to  be  relatively  in- 
significant in  view  of  the  limited  data  base.  (The  authors  emphasized 

that  " . the  sizes  of  the  data  samples  involved  are  small.") 

The  variations  of  the  distribution  for  the  two  aircraft  targets  appears 
to  be  on  the  order  of  20%. 

The  amplitude  vs.  time  histories  were  presented  for  the  repeater 
and  aircraft  over  a total  duration  of  14  minutes.  During  port  ions  of  this 
period  the  repeater  and  aircraft  target  were  present  simultaneously. 

A visual  inspection  of  this  data  did  not  indicate  any  significant  dif- 
ferences in  the  characteristic  of  the  amplitude  variation  of  the  repeater 
and  the  aircraft  echoes.  In  fact,  during  the  major  portion  of  the  time 
when  the  echoes  were  observed  simultaneously, ( which  extended  over  a 
period  of  some  5 minutes)  the  two  sets  of  echoes  appeared  to  fade  with 
similar  periods. 


On  the  basis  of  these  results,  it  would  appear  that  the 
motion  of  the  aircraft  had  little  effect  on  the  fading  characteristics 
observed  by  an  OTH  radar.  This  would  imply  either  that  the  observed 
fading  represents  a temporal  variation  with  relatively  little  spatial 
structure  or  that  the  fading  rate  was  too  rapid  to  be  unambiguously 
represented  by  the  one-second  samples  taken  with  the  radar--so  tnat  the 
observed  fade  rate  is  system,  rather  than  phenomenon,  limited. 


Intermediate  Summary  and  Conclusions 


To  summarize  the  investigation  of  the  possibility  of  using 
fading  characteristics  as  a propagation  mode  discriminant,  the  following 
points  may  be  made. 

1.  The  amplitude  distributions  that  have  been  observed  are 
similar  for  signals  propagated  via  the  E and  F modes. 

2.  Irregularity  drift  measurements  indicate  similar  drift 
(or  wind)  speeds  occur  in  the  two  layers. 

3.  The  mean  irregularity  size  deduced  from  spaced  receiver 
measurements  indicate  that  the  F layer  irregularities  have  a mean 
scale  size  about  twice  that  of  the  F.  layer  but  that  there  is  a sub- 
stantial overlap  in  the  size  distributions  that  were  observed. 

4.  Irregularity -induced  fading  occurs  a substantial  fraction 
of  the  time  (>50%  of  the  time  in  the  Buckland  Park  data)  and  a 
practicable  mode  identification  technique  must  be  operable  in  their 
presence;  this  precludes  complete  dcp.ndence  upon  the  Faraday  rotation 
effect  for  mode  identification. 


i We  may  conclude  from  the  above  that  although  the  fading  pro- 

I 

duced  by  ionospheric  irregularities  on  signals  may  differ  somewhat  in 
scale,  size  and  frequency  of  occurrence,  depending  upon  the  ionospheric 
l layer  supporting  the  propagation  of  the  signal,  the  differences  are 

i 

statistical  rather  than  deterministic.  For  this  reason,  the  fluctua- 
! tion  of  signal ampl i tudes  do  not,  of  themselves,  provide  a conclusive 

indication  of  the  mode  supporting  propagation  of  echoes  from  a given 
target.  Thus,  while  amplitude  fluctuations  provide  a diagnostic  means 
that  may  be  helpful  in  determining  mode  structure  (i.c.,  the  distribu- 
tion of  mode-observable  regions  in  the  f-T  plane)  they  are  of  little 
value  by  themselves  in  determining  the  propagation  mode  of  specific 
detections. 


» 
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5.  ALGORITHM  FOR  MODE -IDENTIFYING  TRACKERS 


This  section  discusses  the  incorporation  of  the  previously 
discussed  mode  identification  criteria  into  an  automatic  tracking 
algorithm.  These  considerations  include  the  provision  of  necessary 
information  regarding  mode  structure  to  the  tracker,  the  information 
to  be  retained  by  the  tracker  concerning  mode  structure,  and  the  logic 
of  the  decision  processes  to  be  applied. 

The  discussion  will  assume  the  use  of  a predictive  filter  for 
tracking.  The  description  of  the  algorithm  and  a justification  for  its 
use  are  presented  in  the  first  section.  The  following  section  dis- 
cusses the  modifications  required  to  make  the  tracker  apply  the  boundar) 
crossing  criteria,  and  the  final  section  discusses  a second  technique 
for  collapsing  multiple  modes  by  Kalman  filtering. 

5.1  Considerations  of  the  Tracking  Filter 

The  problem  of  automatic  detection  and  tracking  of  targets  by 
the  use  of  HF  OTH  radar  is  difficult  and  complex.  Procedures  are 
needed  for  making  detections,  converting  the  coordinates  of  these  de- 
tections (usually  time  delay,  doppler  shift,  and  electrical  phase 
difference)  to  geographical  coordinates  (range,  range  rate,  and  acimuth), 
and  smoothing  the  resulting  tracks.  Figure  5-la  shows  a block  diagram  of 
a straightforward  method  of  carrying  out  these  tasks. 

In  practice,  the  system  described  by  Figure  5-lb  has  the  potential 
to  outperform  that  of  Figure  5-ia.  The  tracker  can  feed  back  information 
to  the  detector,  which  will  aid  the  detector  by  telling  it  where  to 
look  for  detections  during  the  next  time  frame  of  data.  The  coord  mat  •• 
conversion  procedure  can  now  utilize  both  the  raw  and  smoothed  data  in 
making  its  decisions  and  calculations.  The  tracker  will  probably  re- 
quire a preliminary  type  of  coordinate  conversion  on  its  input  data.  \n 
additional  feedback  path  from  the  coordinate  conversion  block  to  the 
tracking  block  may  also  be  required,  although  this  path  is  not  shown  in 
figure  5- lb. 
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(a)  System  Usin»>  No  feedback 


(b)  System  IJsinji  feedback 


I- i "lire  1 . I wo  Possible  System  Configurations 
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The  ultimate  task  is  to  choose  the  type  of  tracking  filter  which 
should  be  utilized  and,  unavoidably  but  to  a lesser  extent,  how  this  filter 
should  interact  with  the  detection  and  coordinate  conversion  routines.  Choos- 
ing a particular  type  of  tracking  filter  from  the  wide  variety  of  available 
filters  is  not  an  easy  task.  Both  aerospace  and  automatic  control  literature 
abound  with  articles  on  Wiener  and  Kalman  filters  (over  200  articles  are  listed 
in  the  bibliography  of  Kailath,  1974).  The  alpha-beta  filter  (Singer  and  Behnke, 
1971)  is  probably  the  next  most  widely  discussed  filter  and,  because  of  its  sim- 
plicity, has  probably  been  implemented  more  often  in  tracking  systems  than  any 
other  filter.  A recent  report  by  Cantrell  (1975)  discusses  adaptive  versions  of 
alpha-beta  and  alpha-beta-gamma  trackers.  Filters  which  achieve  smoothing  by 
fitting  polynomials  to  the  observed  data  have  also  been  proposed  (Suyemoto, 

1973) .  Recently,  nonlinear  adaptive  trackers  (Bershad,  Merryman  and  Sklansky, 

1974)  have  been  investigated;  but  these  results  are  in  the  formative  stage 
and  have  not  been  considered  here. 


Singer  and  Behnke  provide  a detailed  comparison  of  the  Wiener  filter, 
the  alpha-beta  filter,  and  the  Kalman  filter.  The  relative  performance, speed 
and  memory  requirements  of  particular  implementations  of  each  of  these  track- 
ers is  discussed.  The  alpha-beta  filter  assumes  that  the  target  moves  in  a 
straight  line.  Although  maneuver  detection  devices  can  be  used  to  improve 
this  filter's  capability  to  track  maneuvering  targets,  these  devices  are 
probably  not  applicable  to  real-time  HF  radar  systems  since  the  detection 
of  maneuvers  at  a low  signal-to-noise  ratio  and  a sampling  interval  of  one 
second  or  more,  is  itself  a very  difficult  problem.  Thus,  one  criterion 
for  an  HF  radar  tracking  algorithm  should  be  the  ability  to  follow  maneuver- 
ing targets  without  the  need  for  a separate  maneuver-detection  algorithm. 

Two  filters  which  should  be  considered  for  the  application  arc  the 
Kalman  filter  and  the  Wiener  filter.  The  Wiener  filter  is  actually  a sim- 
plified type  of  Kalman  filter  in  which  the  gain  vector  (which  determines 
how  fast  the  tracked  coordinates  can  change)  does  not  vary  with  time  and  is, 
in  fact,  equal  to  the  steady  state  gain  vector  of  the  corresponding  Kalman 
filter.  In  a Kalman  filter,  the  calculation  of  the  gain  vector  requires  a 
considerable  portion  of  the  total  computation  time.  In  the  Wiener  filter, 
the  gain  vector  is  calculated  off-line  and  is  used  repeatedly  thereafter. 
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In  cases  for  which  the  Kalman  filter  rapidly  reaches  the  steady  state,  the 
performance  of  the  Wiener  filter  will  closely  approach  that  of  the  Kalman 
filter  while  requiring  about  one-fifth  the  computation  time  (Singer  and  Behnke, 
19711. 


The  Kalman  filter  has  an  additional  advantage  over  the  Wiener  filter. 
In  the  course  of  calculating  the  variable  gain  vector,  the  Kalman  filter 
computes  an  estimate  of  the  covariance  matrix  of  the  filter's  prediction 
error.  This  matrix  provides  useful  information  about  the  probability  dis- 
trioution  of  the  prediction  error.  If  the  prediction  errors  can  be  assumed 
to  be  zero-mean  random  variables  with  a jointly  Gaussian  distribution,  the 
covariance  matrix  completely  specifies  their  joint  distribution  function. 

This  distribution  function  can  be  fed  back  to  the  detector  (along  with  the 
predicted  values)  and  used  to  make  improved  detections  on  the  next  frame  of 
data  (Petting.  1974) . 

For  our  present  purpose,  we  will  assume  the  more  general  Kalman  filter 
will  be  used  in  the  HF  radar  tracking  scheme  shown  in  Figure  5- lb.  However, 
this  is  done  with  the  reservation  that  should  future  testing  indicate  that 
the  simpler  Wiener  filter  can  produce  results  comparable  to  those  of  the 
Kalman  filter,  it  should  be  used  in  order  to  effect  significant  savings  in 
computer  time  and  program  complexity.  A Kalman  filter  tracking  system  has 
been  tested  and  evaluated  by  Oetting.  Insofar  as  possible,  such  testing 
should  be  made  as  specific  to  the  situation  as  is  possible  to  insure  a valid 
comparison  of  filter  performance. 
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The  following  material  in  this  section  is  taken  from  Petting 
(1974)  to  lay  the  groundwork  for  later  discussions. 

In  order  to  understand  the  origin  of  the  Kalman  filtering 
equations,  we  consider  a one-dimensional  example  of  target  motion,  in 
which  only  xr(t),  the  target’s  range,  is  of  interest.  We  can  write  a 
Taylor  series  expansion  for  xr(t)  in  terms  of  its  derivatives  as: 

xr(fT)  = xr(t)  + Txr(t)  + *sT2kr(t)  ♦ u(t)  (1) 

where  T is  an  interval  of  time  which  is  small  compared  to  the  reciprocal 
of  the  bandwidth  of  xr(t),  and  u(t)  is  the  error  term,  which  results  from 
deviations  in  the  target's  acceleration  ( xrii  0).  Thus,  u(t)  is  usually 
referred  to  as  the  maneuver  noise.  The  Kalman  filter  accepts  imperfect 
measurements  of  xf(t)  and  xr(t)  and  produces  smoothed  estimates  of  xr(t) 
*r(t),  and  k"r(t).  These  measurements  are  presumed  to  be  corrupted  by 
additive  noise,  as  follows: 

yr(t)  = xr(t)  ♦ vr(t)  (2) 

y*(t)  = xrct)  + vt(t) 

where  yr  and  are  the  measured  range  and  range  rate,  and  vr  and  v^ 
are  the  measurement  errors. 

In  order  to  apply  Kalman  filter  theory  to  the  problem  at  hand, 
it  is  necessary  to  write  the  discrete  versions  of  equations  (1)  and  (2): 

xr(k=l)  = xr(k)  ♦ Txr(k)  * hl2XTW  + u(k)  (3) 

yr(k)  = xr(k)  + vr(k) 

= xr(k)  ♦ Vf(k)  (4) 

where  k « 0,  1 , 2,  . . . is  the  index  for  the  sampled  values  of  each  time 
function,  and  T is  now  the  sampling  interval.  For  the  case  where  u(k), 

\ ( k ) , and  v^(k)  arc  sampled  Gaussian  white  noise  processes,  it  can  he 


shown  that,  given  the  measured  quantities  yy  and  y^,  the  Kalman  filter 
calculates  estimates  for  xr,  x^.,  and  x„  (denoted  by  xr,  x^,  and  x.r.) 
which  are  optimal  in  the  sense  of  least  mean  square  error,  where  the 
range  error  er  = yf  - xf  with  similar  expressions  for  range  rate  and 
acceleration  errors. 

While  the  assumption  that  the  measurement  noises  vr  and  v . 
are  white  and  Gaussian  is  a good  assumption  in  most  cases  of  interest  to 
us,  the  maneuver  noise  u will,  in  general,  be  neither  white  nor  Gaussian 
(Kailath,  1974).  This  situation  can  be  handled  in  two  ways  --  we  can 
accept  the  suboptimality  of  the  Kalman  filter  for  non-white  maneuvers,  or 
we  can  whiten  the  maneuver  noise  so  that  optimal  filtering  can  be  achieved. 
Because  of  the  considerable  added  computational  complexity  incurred  by 
the  whitening  procedure,  it  may  be  desirable  to  accept  the  suboptimality 
of  the  Kalman  filter  rather  than  implementing  the  so-called  extended  Kalman 
filter.  Of  course,  in  the  case  where  the  target  exhibits  constant  radial 
velocity  or  constant  radial  acceleration,  the  simple  Kalman  filter  will 
achieve  optimality. 

5.2  Simplified  Two-Dimensional  Kalman  Filter 

We  will  now  discuss  the  two-dimensional  (range  and  azimuth) 

Kalman  filter  which  we  propose  for  use  in  an  HF  radar  system  capable 
cf  measuring  range,  doppler,  and  azimuth.  It  will  be  convenient  to 
use  matrix  notation.  Since  the  independently  measureable  quantities  arc 
range,  range  rate,  and  azimuth,  the  state  vector  x(k)  is  given  by: 

x(k)  = [xr(k),  x^k),  xy  (k),  xQ(k),  x^(k)]r  (5) 

where  xn  and  x.  are  the  azimuth  and  its  time  derivative,  and  the  t 

e e 

indicates  transpose.  The  equations  of  state  are: 

x(k+l)  = <5>x(k)  + Gu(k)  (6) 

y(k)  = Hx(k)  + v(k) 

where  u(k)  = [ur(k)  , u^k)]1  (7) 
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is  the  maneuver  noise  vector 


y£k)  = [yr(k),  yf.Cic),ye(ky] t 


(8) 


is  the  measurement  vector. 


v(k)  » [v  (k)  , v^(k)  , VgCkJ]1 


(9) 


is  the  measurement  noise  vector,  and 


<J>  = 


1 

0 

0 

0 

0 


T hi*-  0 0 

1 T 0 0 

0 10  0 
0 0 1 T 

0 0 0 1 


(10) 


is  the  transition  matrix.  The  matrices  G and  H are  given  by: 


G 


0 0 
0 0 


and  H 


1 0 
0 1 
0 0 


0 0 0 
0 0 0 
0 1 0 


(11) 


At  the  K-th  sampling  instant,  the  Kalman  filter  accepts  the  vector  y(k) 
and  recursively  calculates  four  quantities:  an  estimate  of  the  covariance 

matrix  of  the  prediction  error  at  time  k based  on  all  samples  processed 
prior  to  time  k (denoted  by  P (k | k- 1 ) ; a filtered  estimate  of  x(k)  based 
on  all  simples  processed  thusfar  (denoted  by  £(k|k)  );  an  estimate  of 
the  covariance  matrix  of  the  estimation  error  based  on  all  samples  which 
have  been  processed  (denoted  by  P(k|k)  );  and  a prediction  for  x(k+l) 
based  on  samples  processed  through  time  k (denoted  £(k+l)  ). 


The  expressions  which  must  be  evaluated  (Kailath,  1974)  are: 


s 

I 

I 

i 


P(kjk-l)  - ^(k-llk-l)^  + GQG* 

x (k 1 k)  = x(k|k-l)  + P(k|k-l)HtW-1(k)  [y(k)  - Hx(k|k-lj] 
P(k|k)  = P (k | k- 1)  - P (k | k- 1 ) 1 (k)HP(k|k-l) 

5c(k+l  |k)  - «(k  Jk) 


(12) 

(13) 

(14) 

(15) 


In  equations  (12) 


(15),  the  previously  undefined  quantities  are: 


(16) 


W(k)  - H PCkjk-l)^  + R 


(17) 


and 


(18) 


2 ° 

where  if)  “ and  ifu*  are  the  variances  of  the  (presumed  Gaussian  and  white) 
r ® 2 2 

maneuver  noises:  and  o , a.  , and  oQ  are  the  variances  of  the  Gaussian 

r r 0 

white  measurement  noises.  In  order  that  the  Kalman  filter  be  optimal, 
these  variances  must  be  known  beforehand;  in  practice,  reasonable  values 
are  chosen  by  experimenting  with  different  combinations  of  values  until 
good  results  are  obtained  in  a wide  variety  of  experimental  scenarios. 
Although  a dependence  on  k can  be  introduced  into  these  parameters, 
they  are  indicated  as  constants  in  the  above  equations,  since  we  lack 
prior  knowledge  of  their  possible  time  variations. 
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5.3  Modification  of  the  Tracker  Algorithm  for  Mode 


Identification 

From  the  preceding  discussion,  we  can  identify  the  following 
information  that  can  be  available  about  each  track  maintained  by  the 
automatic  tracker. 

a.  A set  of  unsmoothed  positions  in  the  radar  observables 

space  at  which  echoes  used  to  develop  the  track  were 
detected.  Each  observation  comprises  the  following 
information:  radar  delay,  doppler  shift,  azimuth, 

fine  azimuth,  signal  level  (correlator  output), 
frequency  of  observation  and  time  of  the  observation. 

b.  A set  of  smoothed  estimates  of  position,  each  estimate 
comprising  the  same  quantities  as  those  given  above. 

c.  A set  of  smoothed  predictions  of  future  positions 
in  the  mode  observable  region  at  which  echoes  are 
anticipated,  each  prediction  comprising  the  same 
quantities  as  those  given  above. 

d.  Estimates  of  the  covariance  matrix  of  the  prediction 
errors  and  of  that  of  the  estimation  errors. 

Normally,  an  operating  tracker  would  not  maintain  all  of  this 
information.  Instead,  it  would  drop  the  raw  observations  y(k)  (see 
Figure  5-2)  and  retain  for  display  only  the  set  of  smoothed  estimates 
of  the  targets  position  x (k|k)  (all  k up  to  present),  the  most  recent 
prediction  of  position  x (k+l|k),  and  the  most  recent  estimates  of  the 
covariance  matrix  of  the  estimation  error  P(k|k). 

We  now  consider  the  application  of  the  boundary  crossing 
criteria  and  assume  an  m of  n decision  logic  is  to  be  used.  At  first 
we  will  consider  the  information  and  algorithm  required  when  the  mode 
identifying  tracker  is  to  be  used  with  a frequency  hopping  radar. 

Although  this  mode  identifying  technique  is  conceptually  quite 
simple,  its  implementation  becomes  somewhat  complicated  by  the  need  to 
defer  decisions  until  a sufficient  number  of  observations  has  been 
accumulated.  The  following  four  entities  are  required  within  the  data 
set  maintained  by  the  tracker. 


x (k|k-l)  P(k-1 |k-l) 


Begin 

(k+l)th 

cycle 


Figure  5-2.  Flow  of  Computation  for  Kalman  Filtering. 

Numbers  in  parentheses  are  equation  numbers 
in  text.  All  other  quantities  than  those 
shown  here  are  treated  as  constants. 
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Item  1 


The  time  delay  limits  of  the  various  mode  observable  regions 
(MOR’s)  must  be  recomputed  for  each  frequency.  The  limits  divide 
the  observed  range  of  radar  delay  into  a number  of  segments.  With  regard 
to  a particular  propagation  mode,  each  segment  must  be  classified  as 
one  of  three  hinds:  i.e.,  as  a segment  preceding  the  MOR,  within  the 
MOR,  and  following  the  MOR.  If  the  MOR  is  softly  bounded,  there  will 
also  be  an  interval  of  uncertainty  which  for  the  purpose  of  the  mode 
identification  logic  is  regarded  as  yielding  no  mode  identifying  informa- 
tion. Note  that  each  segment  must  be  classified  as  preceding,  within, 
or  following  the  MOR  of  each  possible  propagation  mode,  and  that  there 
will  in  general  be  several  segments  falling  in  each  category  for  any 
given  mode.  See  Figure  5-3.  One  set  of  this  information,  updated  for 
changed  frequency  on  a dwell-by-dwell  basis,  is  required. 

Item  2 

For  each  segment  of  the  radar  delay,  there  must  be  an  identifica- 
tion of  its  relationship  to  each  MOR.  Two  bits  are  required  to  indicate 
the  relationship  of  a single  segment  to  a single  MOR.  This  information, 
like  Item  1,  can  be  applied  to  all  tracks  and  so  need  only  be  computed 
once  for  each  dwell.  Furthermore,  unless  the  frequency  of  the  radar 
changes  so  much  that  a boundary  for  one  MOR  crosses  that  of  another  MOR 
within  the  interval  of  frequency  change,  the  same  information  can  be 
used  from  one  frequency  dwell  to  the  next.  For  the  present,  we  will 
assume  that  the  further  complication  that  arises  when  MOR  boundaries  cross 
can  be  avoided  and  that  the  table  relating  the  segments  of  radar  delay 
to  the  various  MOR's  can  be  computed  once  for  all  frequencies  of  operation. 

Item  3 

For  each  track,  and  for  each  segment,  there  must  be  maintained 
a record  of  the  successes  and  failures  of  attempts  to  detect  the  echoes 
which  a track  comprises.  For  the  2 of  5 decision  scheme  being  discussed 
here,  these  records  must  include  one  bit  to  indicate  prior  detections  and 
two  bits  to  keep  count  of  prior  "failures-to-detect ."  See  Figure  5-3. 
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Item  4 
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A summary  record  of  decisions  is  made  on  the  basis  of  the  above 
information.  Such  a record  will  be  carried  for  each  of  the  segments 
for  each  track.  Each  record  will  comprise  two  bits:  one  bit  indicates 
a decision  that  echoes  have  been  detected,  and  the  other,  that  echoes 
are  absent  within  the  associated  segment.  Thus,  we  have 


0 0 indicate  insufficient  observations  for  decision 

1 0 M detections  have  occurred 

01  n attempts  resulted  in  fewer  than  n detections 
1 1 not  used. 


to  give  an  indication  of  mode  identity  when  a decision  has  successfully 
been  made  and  of  ambiguity  prior  to  identification. 


Note  that  different  Items  2,  3,  and  4 are  required  for  each 
track.  Item  1 is  recomputed  on  a dwell-by-dwell  basis  and  a single  set 
of  these  limits  applies  to  all  tracks.  Furthermore,  Items  2 and  3 needto 
be  maintained  for  a track  only  until  such  time  that  a mode  identifica- 
tion has  been  made,  after  which  Item  4 suffices. 

Figure  5-3  shows  the  required  information  for  mode  identification 
using  a 2-of-5  detection  scheme.  This  will  be  used  as  an  example  for 
the  subsequent  discussion  of  tne  decision  algorithm.  The  need  for 
the  various  entities  will  become  apparent  in  the  course  of  the  discussion 

The  Kalman  filter  algorithm  yields,  for  each  dwell,  a predicted 
position  for  each  of  the  tracks  being  maintained.  This  position  would 
be  passed  to  the  automatic  detection  processor  which  would  search  the 
incoming  radar  data  for  an  indication  of  the  extension  of  each  track. 
Because  the  present  considerations  of  the  boundary  crossing  criteria 
for  mode  identification  do  not  involve  doppler  information  the 
discussion  will  be  limited  to  the  radar  delay  component  of  the  problem. 

The  automatic  tracker  searches  a neighborhood  of  the  predicted  position. 
For  purpose  of  example,  we  assume  that  this  neighborhood  includes  three 


range  cells  beyond,  and  three  before,  the  predicted  position,  as  well  as 
the  predicted  cell  itself.  If  a signal  exceeding  the  detection  threshold  is 
present  in  any  cell,  an  additional  detection  is  recorded  and  used  to  update 
the  track's  smoothed  estimate,  prediction  and  covariance  matrix.  Note 
that  different  detection  thresholds  can  be  applied  to  different  cells, 
depending  upon  their  position  within  the  neighborhood. 

If  a detection  is  made,  its  position  may  be  compared  with  the 
entries  of  the  limit  table  (Item  1)  to  determine  in  which  radar  delay 
segment  the  detection  has  occurred.  This  can  be  done  by  comparing  the 
ladar  delay  of  the  detection  with  the  radar  delays  listed  in  Item  2 in 
order  of  increasing  delay.  When  a delay  which  exceeds  that  of  the  detec- 
tion is  found,  the  detection  can  be  assigned  to  the  corresponding  inter- 


In  the  event  that  no  detection  is  made  in  the  course  of  search- 
ing the  neighborhood  around  the  predicted  position  of  the  echo,  it 
is  again  necessary  to  assign  a segment  to  the  search  so  that  the  count 
of  the  number  of  attempted  detections  can  be  augmented.  For  this  purpose 
the  predicted  position  is  compared  with  the  limits  of  the  segments.  If 
the  search  neighborhood  is  rather  large,  it  may  be  desirable  to  au,  lent 
the  court  of  attempted  detections  for  all  segments  that  have  more  than 
some  minimal  number  of  cells  within  the  search  neighborhood;  however, 
here  we  will  increment  the  count  only  for  the  region  containing  the 
predicted  position.  Table  5-1  specifies  v...  algorithm  for  performing 
the  above  functions. 

As  the  segments  in  which  a new  detection  or  a failed  attempt  to 
detect  are  identified,  the  information  about  the  region  is  updated.  In 
Table  5-1,  the  index  J identifies  the  segment.  The  modification  of  the 
record  (Items  3 and  4)  mentioned  in  Step  6 of  Table  5-1  requires  only  an 
incrementing  of  appropriate  counters  and  a test  to  see  whether  a decision 
based  on  tr.o  m of  n criterion  can  be  reached.  For  the  purpose  of  this 
discussion,  we  will  refer  to  the  record  for  the  Jth  segment  of  radar  delay 
that  is  maintained  for  detections  which  are  part  of  the  Ith  track  by  the 
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Table  5-1.  Algorithm  for  Determining  Range  Segment  of  a Detection  or 
of  a Search  without  Detection 


N - number  of  possible  modes 


(1)  If  a detection  has  occurred,  set 


D - 1 

R = t (radar  delay  of  detection) 


If  a complete  search  of  the  neighborhood  indicates  no  detection  set 
D = 0 

R = P (predicted  radar  delay) 


W 


K 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


Put  counter  J = 1. 

Compare  R with  limit  T(J) 
Increment  J.  If  J > J 


If  R < t(J)  jump  to  (6). 
, go  to  next  track. 


If  J < J , 
— max 


maximum 
,,  jump  to  (3). 

Modify  record  for  segment  (J)  for  track  under  consideration. 
Go  to  next  track. 


e 


symbol  C(I,J).  Also,  we  make  use  of  the  indicator  D which  is  set  to  1 
if  a detection  has  been  made,  and  0 if  the  search  of  the  neighborhood 
around  the  predicted  position  of  the  detection  indicates  no  detection. 
Figure  5- 3a shows  a flow  diagram  bf  a subroutine  for  performing  the  record 
updating  function. 


If  updating  the  record  does  not  lead  to  a conclusive  result, 
a decision  is  deferred  until  the  results  of  more  observations  are  avail- 
able. If,  on  the  other  hand,  the  new  information  permits  a decision  to 
be  made,  that  the  m detections  in  n attempts  either  has  been  or  cannot 
be  achieved,  a new  attempt  can  be  made  to  identify  the  propagation  mode. 

Identification  of  the  propagation  mode  proceeds  in  the  follow- 
ing stages.  First,  wherever  possible,  modes  are  ruled  out,  either 
because  a detection  (m  of  n)  has  occurred  in  a segment  of  the  radar 
delay  range  that  is  not  within  the  MOR  of  the  eliminated  mode,  or  else 
because  of  a failure  to  detect  (in  of  n)  the  track  within  a segment 
within  the  MOR.  After  as  many  modes  as  possible  are  ruled  out  in  this 
manner,  an  attempt  is  made  to  determine  a sole  surviving  possibility. 

Note  that  the  segmentation  of  the  radar  delay  range  causes  the  criterion 
of  termination  of  a trace  at  a mode  boundary  to  be  applied  simultaneously 
with  that  involving  occurrence  of  the  mode  in  a forbidden  region. 

Figure  5-4  is  a flow  diagram  indicating  how  the  data  regarding 
segmental  detections  and  non- detect ions  can  be  abstracted  for  use 
in  the  decision  process.  For  each  track,  the  records  for  each  segment 
are  examined  to  determine  whether  a decision  has  been  made.  If  so, 
the  segment  map  (Figure  5-3,  Item  2)  M(J,K)  is  examined  tc  determine  the 
implications  of  a decision.  The  inference  can  be  that  the  Kth  mode 
is  ruled  out  or  that  the  detection  (or  failure  to  detect)  is  consistent 
with  the  hypothesis  that  the  Kth  mode  supported  the  propagation  of  the 
echo.  The  decision  that  any  given  mode  (e.g.,  the  Kth)  can  be  ruled 
out  is  recorded  by  setting  S(K)=1.  When  the  implications  of  all  seg- 
mental detections  have  been  drawn  regarding  all  propagation  modes,  the 
algorithm  goes  to  the  final  stage.  This  consists  merely  of  going 
through  the  components  of  S(K)  to  determine  whether  any  of  them  are 
zero.  If  only  one  is  zero,  the  propagation  mode  is  uniquely  identified. 
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Figure  5-3  a.  Subroutine  for  Updating  Record  C(I,J) 
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Abstracting  Process  for  Mode  Identification 
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If  more  than  one  is  zero,  the  identity  of  the  mode  is  still  ambiguous, 
and  further  attempts  at  mode  identification  are  necessary.  If  none  is 
zero,  than  all  modes  have  been  ruled  out.  This  indicates  that  the 
trace  of  the  track  in  the  f-T  plane  is  not  consistent  with  the  mode 
structure  assumed  to  exist.  Aside  from  errors  in  the  assumed  mode 
structure,  such  effects  can  occur  from  errors  in  the  detection  process 
(of  either  kind),  from  targets  not  near  the  earth’s  surface,  or  possibly 
from  ECM  spoof  targets. 

A flow  diagram  for  the  final  stage  of  the  decision  process  is 
shown  in  Figure  S-S. 

5.4  Some  Further  Considerations 

The  preceding  section  discussed  the  algorithms  required  to 
implement  the  mode  identification  criteria  based  on  the  mode  observable 
regions,  as  used  for  a frequency  hopping  radar.  However,  one  problem 
was  neglected- -that  of  indicating  the  case  in  which  the  frequency  changed 
so  much  that  the  boundaries  of  two  MOR’s  cross  within  the  band  of  fre- 
quency- hopping.  Because  the  identity  of  segments  cannot  be  preserved 
across  such  an  intersection  of  boundaries,  a means  is  required  to 
handle  such  cases. 

Figure  5-6  illustrates  an  f-T  plane  in  which  such  a crossing 
of  boundaries  has  occurred.  The  crossing  occurs  at  radar  delay  xc 
and  frequency  f where  the  lagging  boundary  of  the  B MOR  intersects 
the  leading  edge  of  the  A MOR.  Also  included  are  two  mapping  arrays, 
labelled  M(I,J)  and  M*  (I , J) . The  former  of  these  describes  the  relation- 
ship to  the  MOR's  of  the  various  segments  defined  at  frequencies  above 
f , where  the  boundaries  cross  the  latter,  at  frequencies  below  f ^ . 

Note  that  these  two  mapping  arrays  differ  from  one  another. 

Because  the  boundaries  of  MOR's  generally  slope  upward  to  greater 
radar  delay  as  frequency  is  increased,  the  relationship  of  the  segments 
into  which  the  radar  delay  axis  is  divided  at  frequencies  greater  than 
f bear  the  same  relationship  to  the  MOR's  as  do  the  segments  into  which 
the  trace  of  a target  track  in  the  f-T  plane  at  ranges  greater  than  -- 
the  radar  delay  of  the  boundary  intersection.  Consequently,  for  tracks 
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i igure  S-5.  Final  Mode  Identification  Process.  Three  outcomes  arc 
possible:  (1)  ambiguity  remains;  (2)  mode  uniquely 

identified;  (3)  no  mode  survived  testing.  In  the  last 
case,  tests  can  be  applied  again,  or  the  track  can  be 
flagged  as  anomalous  and  be  reported. 
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Figure  5-6.  Diagram  of  the  f-T  Plane  for  a Case  in  which  There  is  a Crossing  of 
MOR  Boundaries.  Different  mapping  arrays  M'(I.J')  and  M(I,.J)  have 
heen  defined  for  use  at  frequencies  below  and  above  the  crossing 
frequency,  f , respectively. 


having  a radar  delay  greater  than  tc,  the  mapping  array  M(I,J)  should 
be  used,  and  for  tracks  with  ranges  less  than  t , the  mapping  array 
M'(I,J)  should  be  used. 

Summarizing,  the  problem  of  boundary  crossings  can  be  readily 
handled  by  defining  different  mapping  arrays  to  be  used  with  tracks 
that  lie  in  different  intervals  of  radar  delay.  In  the  example  shown 
in  Figure  5-6,  the  array  M(I,J)  is  used  for  tracks  with  radar  delays 
exceeding  t_  and  M' (I,J),  for  tracks  with  radar  delays  less  than  t,. 

5.5  Mode  Identification  for  a Radar  Operating  at  Fixed  Frequency 

The  preceding  algorithm  assumes  that  the  radar  hops  from 
frequency  to  frequency  so  that  there  occur  numerous  occasions  when  the 
extension  of  the  trace  into  various  critical  segments  can  be  tested. 

If  these  methods  are  to  be  applied  to  a fixed  frequency  radar,  it  is 
necessary  to  wait  for  the  motion  of  the  radar  target  to  bring  it  to  a 
position  at  which  the  presence  or  absence  of  echoes  critically  indicates 
the  operative  mode  of  propagation.  For  any  given  MOR,  this  happens  only 
twice- -at  the  time  the  target  enters  the  region  the  mode  covers  and 
when  it  exits. 

Thus,  for  fixed  frequency  operation,  the  mode  identification 
algorithm  takes  a very  simple  form:  Regions  immediately  beyond  and 

immediately  ahead  of  each  mode's  MOR  must  be  continuously  monitored 
for  the  presence  of  tracks.  Likewise,  a similar  surveillance  must  be 
placed  on  the  portions  of  the  MOR  just  within  each  of  the  two  boundaries. 
If  a target  is  detected  outside  of  the  MOR,  the  corresponding  mode  is 
dismissed  as  a possibility.  If  it  is  not  detected  outside  the  MOR  but 
is  detected  within  the  MOR,  then  the  propagation  mode  must  be  that 
corresponding  to  the  MOR  being  considered. 

In  many  cases  several  MOR's  have  a common  boundary,  i.e.,  when 
modes  involved  are  various  combinations  of  high  and  low  ray  modes. 

In  these  cases,  there  is  only  one  critical  boundary  which,  regrettably, 
is  a soft  boundary.  Consequently,  the  ability  to  identify  modes  in  these 
cases  hinge  entirely  on  the  radar's  ability  to  definitively  determine 
the  presence  or  absence  of  echoes  in  a single  critical  region. 
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5.6 


Modifications  for  Soft  Boundaries 


1 


If  a soft  boundary  occurs,  a region  must  be  delineated  by  a 
boundary  outside  of  which  the  presence  of  echoes  can  be  detected  with 
the  required  certainty  and  another  outside  of  which  the  absence  of 
echoes  can  be  established  with  the  required  certainty.  The  statistics 
of  these  determinations  is  discussed  elsewhere,  as  is  the  choice  of  the 
contours  of  constant  SNR  at  which  these  boundaries  must  be  placed. 

i 

The  newly  defined  region  is  one  in  which  the  presence  or  | 

absence  of  echoes  is  not  a definitive  indication  of  mode  identity.  \ 

These  regions  can  be  used  to  segment  the  radar  delay  axis  in  the  same 
manner  as  the  MOR's  were  so  used.  The  presence  of  such  segments  will 
be  flagged  by  a zero  in  the  appropriate  part  of  the  mapping  array  M(I,J). 

The  presence  of  the  zero  indicates  that  detections  in  corresponding 
segment  convey  no  information  regarding  the  corresponding  mode.  Figure  5-7 
indicates,  respectively,  the  relationship  between  the  distribution  of 
MOR’s  in  the  f-t  plane  in  such  cases  and  the  resulting  mapping  array. 

« 

i 

) 

| 

i 
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Figure  5-7.  The  Soft  Boundary  at  the  Lagging  lidge  of  the  A MOR 

Leans  to  the  Placement  of  a 0 in  the  Segment  4 Position 
of  the  Manning  Array  M 


The  Mode  Collapsing  Process 


When  it  is  possible  to  keep  track  of  the  relative  distance 
separating  echoes  that  come  from  the  same  target,  but  via  different  modes, 
the  possibility  exists  for  the  tracker  to  coniine  the  various  echoes.  Such 
a procedure  will  be  discussed  in  this  section.  Because  the  method  seems 
to  be  particularly  attractive  for  working  with  high-ray/low- ray  combina- 
tions of  modes,  we  will  concentrate  or.  that  case. 

At  the  skip  distance  the  difference  in  radar  delays  between 
the  high,  low,  and  combination  mode  radar  echoes  is  zero.  As  distance 
increases  beyond  the  skip  distance,  the  relative  delay  between  the 
modes  increases  non-lincarly,  but  monotonically.  Let  us  suppose  that  the 
separation  in  radar  delay  between  the  low-ray-mode  propagated  echo  and  the 
mixed -mode -propagated  echo  is  A (t)  where  r is  the  radar  delay  of  the  low-ray- 
mode echo.  Because  of  the  symmetries  of  the  problem,  the  high -ray  mode 
will  lag  the  mixed-mode  also  by  A(t)  and  the  low-mode  by  2 A (x) . We 
also  suppose  that  the  monotonic  function  A(x)  can  be  expanded  as  a power 
series 


A(x)  = ar 


+ Oj  T + a2 


The  Oj's  of  the  above  expression  are  functions  of  time  because  the  relative 
delay  A(x)  will  vary  in  its  form  as  the  state  of  the  ionosphere  changes. 

The  tracker  will  be  made  tolerant  to  errors  in  the  prediction 
of  range  by  designing  it  to  search  out  a finite  area  of  the  radar  observ- 
ables space,  as  we  discussed  in  earlier  section.  Consequently,  the  power 
series  in  t can  be  truncated,  so  long  as  the  resulting  error  in  A(x)  is 
smaller  than  the  tracker's  search  area.  If  a reasonably  accurate  set 
of  estimates  of  the  ct's  can  be  maintained,  the  power  series  approximation 
to  A(x)  can  be  used  to  predict  the  radar  delay  of  an  echo  propagated  by  a 
later  mode,  given  that  it  was  detected  at  a radar  delay  x. 

The  ability  to  do  this  permits  the  effect  of  multiple  propagation 
modes  to  be  accounted  for.  Two  flow  diagrams  of  variations  on  the  Kalman 
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filter  tracker  with  provision  for  multiple  modes  are  shown  in  Figures  5-8 
and  5-9.  In  Figure  5-8  the  Kalman  filter  used  as  a tracker  is  augmented 
with  a second  estimator-predictor  loop.  In  this  second  loop,  the  ex- 
pansion coefficients  a are  estimated.  In  this  scheme  the  detections 
made  by  the  later  modes  do  not  contribute  to  the  maintenance  of  the 
track.  Rather,  it  is  assumed  that  the  target  can  readily  be  tracked  by 
the  first  mode.  The  predicted  target  state  vectors  X (k+l|k)  are  then 
modified  to  yield  predicted  state  vectors  for  the  later  mode  echoes, 
e.g.,  x'  (k+l|k),  making  use  of  a correction  in  radar  delay  calculated 
from  the  estimates  of  the  a's,  the  coefficients  of  the  power  series  ex- 
pansion of  the  mode  separation  A(t).  The  echoes  are  then  sought  in  the 
neightborhood  of  X (k+l|k)  and  deviations  of  the  observations  y'(k+l)  are 
then  fed  back  to  the  estimator,  where  the  estimates  of  the  coefficients 
a are  updated.  The  time  constants  for  the  estimates  of  a must  be  chosen 
to  be  relatively  long  compared  to  those  used  in  the  estimates  of  the  target 
state  vectors  x(k|k),  so  that  target  prediction  errors  do  not  significantly 
affect  the  estimates  of  the  a's. 

Figure  5-9  diagrams  a second  method,  in  whic'..  detections  made  via 
the  later  modes  are  converted  to  the  range  of  the  first  mode  where  they 
can  be  included  as  data  points  for  estimating  the  track.  However,  in 
order  to  improve  the  estimates  of  the  a’s  the  observations  are  separated 
according  to  mode  and  the  a's  are  corrected  to  remove  any  observed 
bias. 

This  schem,  provides  a greater  number  of  data  points  with  which 
the  tracker  can  work  but  does  not  improve  the  ability  to  estimate  the 
a's.  These  are  based  upon  a division  of  the  available  data  points  into 
three  groups,  according  to  propagation  mode.  However  if  3S  before  the 
estimates  of  the  a's  are  based  on  a longer  history  of  observations  than  are 
those  of  the  tracks,  the  time  averaging  can  more  than  make  up  for  the  ef- 
fect of  the  division  of  the  data  points  into  groups  according  to  propagation 
mode . 
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Figure  5-8.  Modification  of  the  Tracker  Algorithm 
extract  Multiple  Mode  effects. 
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A further  note  regarding  the  operation  of  the  tracker  in  Fig- 
ure 5-9  concerns  the  relative  signal  strengths  of  the  echoes  propagated 
via  the  different  inodes.  Usually,  one  of  the  modes  will  provide  , ss 
path  loss  than  does  the  other,  with  the  result  that  there  may  he  sub- 
stantial differences  in  the  signal-to-noise  ratio  of  echoes  received  by 
one  mode  and  those  of  another.  Consequently,  the  probability  of  de- 
tection and  false  alarm  rate  cannot  be  kept  the  same  for  the  detection 
of  echoes  propagated  by  the  different  modes.  This  difference  would  have 
to  be  accounted  for  in  the  tracker  so  that,  for  example,  the  tracker  would 
not  drop  track  on  the  basis  of  having  failed  to  detect  the  echo  as 
propagated  by  a particularly  weak  mode. 


The  above  discussion  has  concentrated  on  the  radar  delay  of  the 
observed  targets.  However,  a similar  set  of  conditions  apply  to  the 
doppler  shift  as  well.  Just  as  was  true  for  the  radar  delay,  there  will 
be  a difference,  from  one  mode  to  another,  of  the  doppler  shift  observed 
on  the  echoes.  Again,  these  can  be  expanded  as  a function  of  radar  delay 
and  used  to  interrelate  the  state  vectors  describing  the  echo  as  observed 
by  the  different  modes. 
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6.  CONCLUSr.S 


The  existence  of  multiple  propagation  modes  can  have  two  kinds 
of  effects  on  the  data  developed  and  displayed  by  an  OTH  radar.  The 
first  effect  is  to  introduce  an  ambiguity  in  the  interpretation  of 
tracks.  Correct  identification  of  the  propagation  mode  is  necessary 
before  the  observed  radar  delay  of  the  echo  can  be  interpreted  in  terms 
of  a distance  from  the  radar.  Misidentification  can  result  in  an 
error  amounting  to  several  tens  to  hundreds  of  km  in  the  position 
assigned  to  an  observed  track.  Because  echoes  propagated  by  different 
modes  can  occur  at  a common  range,  the  risk  of  making  such  an  error 
exists. 

The  second  effect  is  to  cause  a single  aircraft,  or  other  target, 
to  produce  multiple  echoes  at  the  radar,  each  having  propagated  by  a 
different  propagation  mode  or  combination  of  such  modes,  and  consequently, 
each  arriving  with  a different  radar  delay.  If  the  common  origin  of 
these  echoes  is  not  recognized,  the  ■'•adarist  will  overestimate  the 
number  of  targets  he  is  viewing.  Conversely,  if  he  wrongly  associates 
echoes,  assuming  them  to  have  been  caused  by  multimoding  when,  indeed, 
they  have  not  been,  he  will  underestimate  the  number  of  targets. 

Multimoding  between  distinct  propagation  modes  can  occur  during 
the  daytime  when  the  F layer  is  accompanied  by  sunlight  sensitive 
underlying  layers,  namely  the  E and  FI  layers.  It  can  also  occur 
during  summer  nights  when  the  sporadic  E layer  occurs.  However,  con- 
fusion of  E-and  F- layer  propagated  echoes  is  only  a problem  at  the  near 
ranges  reached  by  the  relatively  short  distance  E layer. 

The  risk  of  confusion  also  exists  between  one-hop  and  two-hop 
modes  of  a single  layer.  These  can  usually  be  avoided  if  the  operating 
frequency  of  the  radar  can  be  made  to  exceed  75%  of  the  maximum 
observable  frequency  (MOF)  at  the  range  of  interest. 

Finally,  the  confusion  of  high- ray  and  low- ray  propagated  echoes, 
and  combinations  of  these,  can  also  lead  to  problems.  Because  of 
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defocussing  effects,  the  high  ray  becomes  significantly  attenuated 
at  frequencies  more  than  a couple  Ffiz  below  the  MOF  to  the  range  of 
interest.  Unfortunately,  the  obscuring  effect  of  underlying  layers 
may  preclude  such  a choice  of  operating  frequency. 
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A number  of  tests  can  be  applied  to  determine  the  propagation 
mode  responsible  for  observed  track.  Many  of  these  involve  knowledge 
of  the  maximum  observable  frequency  to  the  range  in  question  and  to  other 
details  of  the  propagation  conditions, such  as  the  maximum  radar  delay 
at  which  echoes  propagated  via  a particular  mode  can  appear. 

To  facilitate  the  discussion  of  these  mode  identifying  tests, 
a number  of  concepts  were  introduced.  Among  these  were  the  mode  obser- 
vable region  (MOR)  for  a given  propagation  mode,  being  the  region  in 
the  frequency-  delay  plane  where  surface  targets  (which  here  includes  low 
aircraft)  can  occur.  The  boundaries  of  the  mode  observable  regions  can 
be  discretely  bounded  or  can  be  the  consequence  of  a gradual  diminution 
of  signal  strength.  The  adjectives  "hard"  and  "soft"  were  applied 
respectively  to  such  boundaries. 

It  was  shown  that  the  detectability  of  echoes  on  either  side 
of  such  boundaries  provided  a means  for  testing  the  identity  of  the 
propagation  mode  responsible  for  any  given  track.  The  statistical 
considerations  required  for  such  a determination  at  hard  and  soft 
boundaries  was  related  to  the  maximum  acceptable  probability  of  mis- 
identifying  the  propagation  mode. 

Various  tracker  algorithms  were  considered,  and  particular 
attention  was  given  to  the  suitability  of  a Kalman  filter  tracker 
for  incorporation  into  a mode  identifying  tracker.  To  be  suitable, 
it  is  necessary  that  the  associated  automatic  detection  scheme  be 
capable  of  producing  a null  decision. 

Algorithms  were  provided  for  applying  the  mode  identifying 
criteria  to  tracks  as  developed  by  a Kalman  filter.  As  an  example, 
the  case  where  the  required  certainty  of  mode  identification  would  be 
achieved  with  a 2-out-of-5  detection  logic  was  considered. 
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Additional  flow  diagrams  are  presented  for  modifying  the  Kalman 
filter  to  perform  a mode- col lapsing  process  whereby  multiple  tracks 
produced  by  a single  target  could  be  associated  by  maintaining  an  es- 
timate of  mode  separation  as  a function  of  radar  delay. 

This  latter  process  has  particular  value  for  the  case  when  the 
multimoding  is  the  result  of  high- ray,  low- ray  and  combination  modes. 
In  such  a case,  the  different  modes  have  one  boundary  in  common,  while 
the  other  boundaries  are  soft,  making  the  application  of  the  boundary- 
crossing  criteria  difficult. 


